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Abstract 

The description of dispersion forces within the framework of macroscopic quantum 
electrodynamics in linear, dispersing and absorbing media combines the benefits of 
approaches based on normal-mode techniques of standard quantum electrodynamics 
and methods based on linear-response theory in a natural way. It renders generally 
valid expressions for both the forces between bodies and the forces on atoms in the 
presence of bodies while showing very clearly the intimate relation between the dif- 
ferent types of dispersion forces. By considering examples, the influence of various 
factors like form, size, electric and magnetic properties, or intervening media on the 
forces is addressed. Since the approach based on macroscopic quantum electrody- 
namics does not only apply to equilibrium systems, it can be used to investigate 
dynamical effects such as the temporal evolution of forces on arbitrarily excited 
atoms. 

Key words: dispersion force, macroscopic quantum electrodynamics, Casimir 
effect, van der Waals force, atom-surface interaction, intermolecular potential, 
atomic polarizability, magneto-electric medium, multilayer structure, spontaneous 
decay, weak atom-field coupling, strong atom-field coupling, Rabi oscillations 
PACS: 12.20.-m, 42.50. Ct, 34.50.Dy, 42.50.Nn 



Contents 

1 Introduction 2 

1.1 Dispersion forces 3 

1.2 Experimental observations 4 



* Corresponding author. 

Email address: s.buhmann@tpi.uni-jena.de (Stefan Yoshi Buhmann). 

URL: www.tpi.uni-jena.de/tpi/qophysics/qo.html (Dirk-Gunnar Welsch). 
1 This work was supported by the Deutsche Forschungsgemeinschaft. 



1.3 Applications 8 

1.4 Theoretical approaches 9 

2 Elements of QED in linearly responding media 19 

2.1 The medium-assisted electromagnetic field 19 

2.2 Atom-field interaction 25 

3 Forces on bodies 32 

3.1 Casimir stress in planar structures 36 

3.2 Macro- and micro-objects 39 

4 Forces on atoms 45 

4.1 Ground-state atoms 45 

4.2 Excited atoms 61 

5 Concluding remarks 75 

Acknowledgements 77 

A Overview over scenarios 78 

B Green tensors 80 

References 81 



1 Introduction 

Dispersion forces originate from the electromagnetic interaction between elec- 
trically neutral objects which do not carry permanent electric and magnetic 
moments. They have been of increasing interest because of their important 
impact on many areas of science. In particular, the extremely miniaturized 
components in nanotechnology can be strongly affected by dispersion forces. 
Recent progress in experimental techniques has led to accurate measurements 
of dispersion forces which have confirmed some of the theoretical predictions 
while posing new questions at the same time. 
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1 . 1 Dispersion forces 



The prediction of dispersion forces is one of the most prominent achievements 
of quantum electrodynamics (QED) where they can be regarded as being a 
consequence of quantum ground-state fluctuations. In order to understand how 
quantum fluctuations lead to dispersion forces, it may be helpful to first recall 
the corresponding classical situation. According to classical electrodynamics, 
electrically neutral, unpolarized material objects will not interact with each 
other, even if they are polarizable. An interaction can only occur if (i) at least 
one of the objects is polarized or (ii) an electromagnetic field is applied to at 
least one of the objects. In the former case the object's polarization will give 
rise to an electromagnetic field which can induce a polarization of the other 
polarizable object (s); in the latter case the applied field induces a polarization 
of the object which in turn gives rise to an electromagnetic field acting on the 
other object(s). Both cases result in polarized objects interacting with each 
other via an electromagnetic field, the interaction and the resulting attractive 
forces between them being a consequence of the departure from the classical 
ground state — unpolarized objects and vanishing electromagnetic field. 

In QED, the state that most closely corresponds to the classical ground state 
is given by the material objects being in their (unpolarized) quantum ground 
states and the electromagnetic field being in its vacuum state, such that both 
the electromagnetic field and the polarization of all objects vanish on the 
quantum average. At first glance, one could hence expect the absence of any 
interaction between the objects. However, the Heisenberg uncertainty principle 
necessarily implies the existence of ground state fluctuations, i.e., both (i) fluc- 
tuating polarizations of the objects and (ii) a fluctuating electromagnetic field 
will always be present. These fluctuations give rise to an interaction between 
the objects — a purely quantum effect which is manifested in the dispersion 
forces acting on them. At finite temperatures, additional thermal fluctuations 
come into play. 

Thus, dispersion forces — also known as Casimir or van der Waals forces 2 — are 
ever-present long-range forces between atoms and/or macroscopic bodies, i.e., 
they exist even if the interacting objects are electrically neutral and do not 
carry electric or magnetic moments. Naturally, dispersion forces have many 
important consequences. On a microscopic level, they influence, e.g., the prop- 
erties of weakly bound molecules [1,2]. A prominent macroscopic signature of 
dispersion forces is the well-known correction to the equation of state of an 
ideal gas, leading to the more general van der Waals equation. 3 But disper- 

2 Often, the term Casimir force is used to denote dispersion forces on a macroscopic 
level whereas dispersion forces on a microscopic level are referred to as van der Waals 
forces. 

3 In fact, it was in this context that the existence of dispersion forces was first 
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sion forces also influence the macroscopic properties of liquids and solids such 
as the anomalies of water [4], the magnetic, thermal and optical properties of 
solid oxygen [5] or the melting behavior of weakly bound crystals [6]. 

The influence of dispersion forces becomes even more pronounced in the pres- 
ence of interfaces between different phases and/or media. Atom-surface dis- 
persion interactions drive the adsorption of inert gas atoms to solid surfaces 
[7-9], influence the wetting properties of liquids on such surfaces [8,10,11] and 
lead to the phenomenon of capillarity [12]. The mutual dispersion attractions 
of colloidal particles suspended in a liquid [13] influence the stability of such 
suspensions [14,15]; unless sufficiently balanced by repulsive forces, they lead 
to a clustering of the particles, commonly known as flocculation [16,17]. 

The above mentioned relevance of dispersion forces to material sciences and 
physical chemistry being rather obvious, it is perhaps more surprising to note 
that they also play a role in astrophysics and biology. Thus, dispersion forces 
initiate the preplanetary dust aggregation leading to the formation of planets 
around a star [18]. Furthermore, they are needed for an understanding of the 
interaction of molecules with cell membranes [19,20] and of cell adhesion driven 
by mutual cell-membrane interactions [19,21]. Recently, dispersion forces have 
been found to be responsible for the remarkable abilities of some gecko [22] 
and spider species [23] to climb smooth, dry surfaces. 

1.2 Experimental observations 

A force between two macroscopic bodies can most easily be measured in a 
quasi-static way where the bodies are brought close together and the force to 
be measured is compensated by a force of known magnitude. In the first use of 
this idea for measuring dispersion forces, the compensating force was provided 
by a Hookean spring and the distance of the bodies was measured by means of 
optical interferometry [24]. In this way, forces between two dielectric [24-28] 
and metal plates [29-32] were investigated. Difficulties in aligning the plates 
were overcome by using alternative setups of a plate interacting with a spher- 
ical lens [28,33-35], two interacting spheres [36], and crossed cylinders [36,37]. 
Magnetic forces generated by electric currents were also used to compensate 
dispersion forces where via feedback, the force values could be inferred from 
current measurements [38-41]. The typical outcome of all these experiments 
was the observation of attractive forces which follow l/z A and 1/z 3 power 
laws for the plate-plate and plate-sphere 4 geometries, respectively, with z 
denoting the object-object separation (see Ref. [42] for a review). However, 
owing to the rather low accuracy — the main limitations being the presence 
of electrostatic forces due to residual charges, the roughness of the samples, 

predicted, for a historical review, see Ref. [3]. 

4 Note that the plate-sphere separation was much smaller than the sphere radius. 
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the lack of accurate position control and the low resolution of the actual force 
measurements — the early results remained controversial. This is best illus- 
trated by the fact that differing power laws [24,25] and even signs [29] were 
found. 

Substantial progress was made by using a torsion-balance scheme [43] where 
very smooth bodies, piezoelectric devices for position control, and capacitive 
force detection, were used to measure the force between a metal sphere and 
plate with high accuracy, thereby confirming an attractive force proportional 
to 1/z 3 . This experiment has been followed by a number of experiments which 
have profited by recent developments in nanotechnology (for an overview, see 
Refs. [44,45]). By attaching a microsphere to the cantilever of an atomic force 
microscope (playing the role of the spring) and monitoring the dispersion- force 
induced bending of the cantilever via deflection of an optical beam, the force 
between the sphere and a nearby surface was measured to obtain very accu- 
rate results for various metals [46-50] and/or dielectrics [51,52], including the 
influence of finite conductivity [46,47,49] as well as surface roughness [48,50]. 
It was further demonstrated that one-dimensional periodic surface corruga- 
tions can lead to a sinusoidally varying tangential force in addition to the 
attractive normal Casimir force [50,53,54]. A similar experimental setup was 
used to measure the force between two dielectric cylinders [55] and spheres 
[56]. Dispersion forces have also been measured by means of a micromachined 
torsion-balance scheme where a small plate suspended by two thin rods ro- 
tates in response to the Casimir force exerted by a nearby microsphere and 
this tilting is monitored by capacitive measurements. The scheme was used 
to study the force between dissimilar metals [57,58] and to demonstrate its 
dependence on the thickness of the interacting objects [59,60]. Changing the 
objects' reflectivity in the visible region by means of hydrogen deposition was 
found to have no observable influence on the force [60,61], indicating that it 
should depend on the frequency-dependent body properties in some integral 
form. 

Dispersion forces can also be measured in a dynamical setup, based on the idea 
that any interaction will affect the relative motion of two objects. In the first 
experimental realization of this idea, a spherical lens was mounted on a loud 
speaker and periodically driven, thereby inducing — by means of the Casimir 
force — a similar motion of a nearby plate mounted on a microphone [62,63]. 
Detection of the amplitude of the induced oscillations led to accurate force 
measurements. This idea has been applied in modern experiments to infer the 
Casimir force from the periodic motion of an object such as a microsphere os- 
cillating at the tip of an atomic force microscope cantilever while interacting 
with a surface [64,65] or an oscillating micromachined torsion balance inter- 
acting with a rigidly mounted sphere [57,58,66,67]. Dynamical measurements 
of this kind can also provide high-precision results for atom-body dispersion 
interactions, which was demonstrated by observing the change of the oscil- 
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latory motion of a single excited ion trapped in a standing electromagnetic 
wave [68] and of cold gases of (ground-state) atoms confined in a magnetic 
trap [69] or an optical lattice [70], induced by their dispersion interaction 
with a nearby surface. In the latter experiment, a temperature dependence of 
dispersion forces was observed. 

Controlled dynamical measurements of dispersion forces on atoms have only 
become feasible recently due to the availability of efficient techniques to cool 
and trap single atoms. In the early experiments, scattering techniques were 
employed which are of course much simpler to implement. Atom-atom inter- 
actions have been observed by scattering a beam of ground-state atoms with 
known narrow velocity profile off a second beam of atoms with thermal veloc- 
ity distribution [71-73] or a stationary target gas [74-77] where an attractive 
1/z 7 force, was found. Scattering techniques have also been employed to probe 
the interactions of atoms with anisotropic molecules [78-80] and even the in- 
teraction of excited atoms with ground-state atoms [81] where in the latter 
case a strong enhancement of the force, was observed. Evidence of atom-body 
forces was first found by observing the deflection of a beam of ground-state 
atoms or molecules passing near the surface of a metal or dielectric cylinder, 
the results suggesting an attractive l/z 4 force [82-85]. 5 In a similar scheme, 
the deflection of atoms passing between two metal plates was monitored by 
observing the atom flux losses due to the sticking of atoms to the plates. In 
doing so, a strong enhancement of the force on excited atoms was observed 
[86]. It was further found [87,88] that the distance dependence of the ground- 
state force changes from a l/z 4 power law for small atom-surface separations 
(non-retarded regime) to the more rapidly decreasing l/z 5 power law as soon 
as the separations exceed the relevant transition wavelengths of the atoms and 
the bodies (retarded regime 6 ). 

It has turned out that introducing a controllable compensating force is also 
useful in the context of atom-scattering experiments; this is the central idea of 
the evanescent- wave mirror: A laser beam is incident on the surface of a dielec- 
tric from the inside at a sufficiently shallow angle, such that total reflection 
leads to an exponentially decaying electric field at its exterior. An atom placed 
in the vicinity of the body will interact with this evanescent field, leading to 
an optical potential. If the laser frequency is larger than the relevant atomic 
transition frequency (blue detuning), then this potential is repulsive; thus cre- 
ating the required compensating force which can be controlled by varying the 
laser frequency and intensity. In this way, dispersion forces on ground-state 
atoms can be measured by monitoring the reflection of the atoms incident on 



In the experiments, the minimum atom-surface separation was so small that to a 
good approximation, the cylinder surfaces can be regarded as planar. 
6 Note that the above mentioned measurements of Casimir forces between macro- 
scopic bodies typically operate in the retarded regime. 
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evanescent-wave mirrors [89-91]. Alternatively, compensating forces can be 
provided by the magnetic fields created by magnetic films, the strength being 
controlled by varying the film thickness [92]. 

Effects due to the wave nature of the atomic motion become relevant for small 
values of the (center-of-mass) momentum such that the atomic de Broglie 
wavelength becomes sufficiently large. In this case quantum reflection of an 
atom from the potential associated with the atom-body force may occur [93]. 
Quantum reflection of ground-state [94-97] and excited atoms [98,99] incident 
on the surface of dielectric bodies was observed in various experiments where a 
detailed measurement of the atom-surface dispersion potential, was achieved 
by recording the reflectivities at different normal velocities. Another prominent 
wave phenomenon that can be exploited for the measurement of atom-body 
potentials is the diffraction of an atomic wave incident on a transmission 
grating forming a periodic array of parallel slits. When passing the slits (which 
may be regarded as small planar cavities), the atomic matter wave acquires 
a phase shift due to the dispersion potential which affects the interference 
pattern forming behind the slits. By comparing the experimental observations 
with theoretical simulations, the interaction of ground-state [100-103] as well 
as excited atoms [104-106] with dielectrics in the non-retarded regime has 
been measured. 

Spectroscopic measurements provide a powerful indirect method for studying 
atom-body dispersion interactions. Here, the fact is exploited that the dis- 
persion potential of an atom can be identified with the position-dependent 
shift of the respective atomic energy level [107]. The resulting shifts of the 
atomic transition frequencies can be observed by spectroscopic means. As the 
shifts are usually much more noticeable for excited levels, this approach yields 
good estimates of the dispersion potentials of atoms in excited energy eigen- 
states. This was demonstrated in experiments measuring dispersion potentials 
of atoms inside planar [108-110] and spherical metal cavities [111,112], near a 
dielectric half space [113,114], and of an ion near a metal plate [115]. In this 
context, selective reflection spectroscopy of atomic gases has proven to be a 
particularly powerful method [116]. It is based on the fact that the reflection of 
a laser beam incident on a gas cell is modified due to the laser-induced polar- 
ization of the gas atoms which in turn is strongly influenced by the dispersion 
interaction of the atoms with the walls of the cell. By comparing measured 
reflectivity spectra with theoretically computed ones, very accurate informa- 
tion on the non-retarded dispersion interaction of atoms with dielectric plates 
[117-122] was obtained, including the potentials of atoms in very short-lived 
excited energy eigenstates which are difficult to study by scattering methods. 
Note that the dispersion interaction with metal plates is much more difficult 
to observe via selective reflection spectroscopy [123]. As a major achievement, 
the method has shown that the dispersion forces on excited atoms can be 
repulsive [124,125]. 
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1.3 Applications 



Taking advantage of the substantially improved sensitivity of dispersion-force 
measurements, comparison of the experimental results with theoretical pre- 
dictions can nowadays even be used to place constraints on other short-scale 
interactions of fundamental interest, such as non-standard gravitational forces 
[44,58,69,126-129]. In addition, dispersion forces have become of increasing im- 
portance in applied science such as nanotechnology and related fields. While 
providing a powerful tool for surface control, e.g., in near-field scanning mi- 
croscopy [130,131], they can also be a disturbing factor whose influence will 
become more and more pronounced with proceeding miniaturization. In par- 
ticular, they can lead to an undesired and permanent sticking of (small) ob- 
jects to surfaces [132-134]. A similarly disturbing effect is observed when atom 
traps are operated near surfaces where dispersion forces can diminish the depth 
of magneto-optical traps, thereby imposing limits upon the near-surface op- 
eration of such traps [135,136]. Traps based on evanescent waves [137-141] 
necessarily operate in the near-surface regime so that dispersion forces au- 
tomatically come into play. The influence of dispersion forces also needs to 
be taken into account when constructing evanescent-wave based elements for 
atom guiding [142-146]. 

Dispersion forces are indispensable in atom optics [147] where mirrors and 
beam splitters for atomic matter waves, have been constructed based on the 
dispersion interactions of atoms with flat surfaces and transmission gratings, 
respectively. Transmission gratings can be used to realize Mach-Zehnder-type 
interferometers for atoms [102]. Flat quantum reflective mirrors provide a fo- 
cussing mechanism when dispersion and gravitational forces are combined with 
in an appropriate way [148]. In addition, by locally enhancing the reflectivity 
of the mirrors via a Fresnel reflection structure [149], reflection holograms for 
atomic matter waves can be realized [150]. The efficiency of atomic mirrors can 
also be enhanced by using evanescent-wave mirrors which can even operate 
quantum-state selectively [151]. As recently reported, the quantum reflection 
of ultracold gases at dielectric surfaces gives rise to interesting phenomena, 
such as the excitation of solitons and vortex structures [152]. 

Further impact on the application of dispersion forces has been made by the 
recent proposal [153] and subsequent fabrication [154] of materials with tai- 
lored magneto-electric properties, also known as metamaterials. 7 Metamate- 
rials displaying simultaneous negative permittivity and permeability in some 
frequency range, allow for the existence of traveling electromagnetic waves 
whose electric-field, magnetic-field and wave vector, form a left-handed triad 
[160] 8 , leading to a number of unusual effects. It is yet an open question 

7 For the current state-of-art of metamaterial fabrication, see, e.g., Refs. [155-159]. 

8 For this reason materials with these properties are commonly referred to as left- 
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whether left-handed properties can lead to interesting phenomena in the con- 
text of dispersion forces and to what extent metamaterials can be exploited 
to tailor the shape and sign of these forces. An interesting behavior of disper- 
sion forces may also occur in conjunction with soft-magnetic alloys, such as 
permalloy or Mumetal [161,162]. After heating and rapid cooling (a process 
called annealing), these materials are in a state of extremely high permeability; 
values of more than 5 x 10 4 have been reported for Mumetal [163]. 



1.4 Theoretical approaches 



As already mentioned, dispersion forces arise from quantum zero-point fluctu- 
ations, namely the fluctuating charge and current distributions of the interact- 
ing objects and the vacuum fluctuations of the (transverse) electromagnetic 
field. If the separation of the objects is smaller than the wavelengths of the 
relevant field fluctuations, then the latter can be disregarded, allowing for a 
simplified treatment of dispersion forces. In this non-retarded regime, disper- 
sion forces are dominated by the Coulomb interaction of fluctuating charge 
distributions. In particular, the Coulomb interaction between two atoms may 
within a leading-order multipole expansion be regarded as the interaction of 
two electric dipoles d and d , 

y_ d-d' -344 m 

4tt£ 2 3 ' { ' 

This approach was first used by London in conjunction with leading-order 
perturbation theory to derive the potential energy of two isotropic ground- 
state atoms to be 

u( ,_c 1 |( jdjfc)|>idW)| 2 

[Z) ° 247^2^ E k + E w - (E + E«) [Z > 

with \k^) and E k (,) denoting the eigenstates and -energies of the unperturbed 
atoms [164]. The London potential implies an attractive force proportional to 
1/z 7 . The idea of deriving dispersion forces from dipole-dipole interactions by 
means of perturbation theory was later applied to three- [165-171], four- [172] 
and A-atom interaction potentials [173-176]. Lennard- Jones showed [177] that 
the interaction of an atom with a perfectly conducting plate can be treated on 
an equal footing, by using the image-charge method. Considering the dipole- 
dipole interaction of the atom with its own image in the plate (instead of a 



handed materials. 



9 



second atom) within first-order perturbation theory, he found a potential 



(0|d 2 |0) 

implying an attractive l/z 4 force. The influence of fluctuating magnetic dipoles 
[178-180], electric quadrupoles [181] and higher multipoles [182,183] as well as 
permanent electric [164,168,184,185] and magnetic dipoles [184] on the atom- 
atom interaction has also been discussed. In this way, it was found that the 
force between a magnetizable and a polarizable atom is repulsive and propor- 
tional to l/z 5 , in contrast to the attractive l/z 7 force between two polarizable 
atoms. Furthermore, studying the interaction of atoms prepared in excited 
energy eigenstates showed that the contributions to the force which arise from 
real, resonant transitions can be attractive or repulsive [168,179,184] (for fur- 
ther reading regarding the atom-atom interaction, refer to Refs. [186-188]). 
Similar extensions have been accomplished regarding the interaction of atoms 
with (planar) bodies. Quadrupole [189] and higher-order multipole atomic mo- 
ments [190-192] were included in the interaction of ground-state atoms with 
perfectly conducting plates and extensions of the image-charge method to the 
interaction of ground-state [193,194] and atoms in excited energy eigenstates 
[195] with planar dielectric bodies were given. 

The method can be further improved by describing the atom and the body on 
an equal footing in terms of their charge densities and expressing the resulting 
interaction potential in terms of electrostatic linear response functions 9 of the 
two systems. This was first demonstrated for a ground-state atom interacting 
with a realistic electric 10 half space exhibiting non-local properties [199]. The 
approach was demonstrated to lead to a finite value of the interaction potential 
in the limit z — > [190,200-203]. 11 For sufficiently large values of z, the 
potential for an atom in front of a half space can be given by an asymptotic 
power series in l/z [199,206-209] 



[a(c<j), dipole polarizability of the atom; e(u), local permittivity of the half 
space] where the leading-order term corresponds to an attractive l/z 4 force 



9 In contrast to the quantities appearing in earlier attempts to treat conductors in 
a more realistic way [196-198], the two response functions are directly accessible to 
measurements. 

10 The term electric is used where no explicit distinction is made between metals 
(conductors) and dielectrics (insulators). Likewise, the notion magneto-electric is 
used to refer to metals or dielectrics possessing non-trivial magnetic properties. 

11 For the dispersion interaction between two dielectric half spaces, this is shown in 
Refs. [204,205]. 
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and coincides with the perfect conductor result (3) in the limit of infinite per- 
mittivity. Next-order corrections are due to the atomic quadrupole polarizabil- 
ity on the one hand [203] and the leading-order non-local dielectric response 
on the other hand [199]. The response-function approach has been used to 
study the interaction of various ground-state objects with half spaces of dif- 
ferent kinds, such as the forces on an ion [210] and a permanently polarized 
atom [211,212] in front of a metal half space, an anisotropic molecule in front 
of an electric half space [213] as well as the interaction of two atoms in front 
of a metal [214,215] and an electric half space [216]. Extensions include the 
interaction of an atom in an excited energy eigenstate with an electric [217] 
and a birefringent dielectric half space [218]; non-perturbative effects [219]; 
effects due to a constant external magnetic field [220]; and the interaction 
of single ground-state atoms/molecules with bodies of various shapes where 
perfectly conducting [221], non-local metallic [222,223] and electric spheres 
[221], non-local metallic [224,225] and electric cylinders [226], perfectly con- 
ducting planar [227] and non-local metallic spherical cavities [228,229], have 
been considered. 

The interaction of two macroscopic bodies B and B' was first treated by 
pairwise summation over the microscopic London potentials (2) between the 
atoms constituting the bodies [230,231], 

^-EZA, (5) 

reBr'eff I 1 1 I 

yielding an attractive 1/z 3 force between two dielectric half spaces [230,232]. 
Though applicable to bodies of various shapes (cf. also Refs. [200,233]), the 
method could only yield approximate results due to the restriction to two- 
atom interactions. By modeling the body atoms by harmonic oscillators, the 
interaction energy of the bodies could be shown to be a sum of all possible 
many-atom interaction potentials [230,232,234]. Applications to the interac- 
tion of two half spaces [175] and two spheres [235] were studied. Microscopic 
calculations of the dispersion interaction between bodies were soon realized to 
be very cumbersome, in particular for more involved geometries. In an alter- 
native approach based on macroscopic electrostatics, the interaction energy 
can be derived from the eigenmodes of the electrostatic Coulomb potential 
which are subject to the boundary conditions imposed by the surfaces of dis- 
continuity [236,237] (for an overview, cf. Ref. [45]). The method was used 
to calculate Casimir forces between electric spheres [237]; electric spherical 
cavities [238]; metal half spaces exhibiting non-local properties [208]; rough 
electric half spaces [239,240]; and electrolytic half spaces separated by a di- 
electric [241]. 

Even though electrostatic methods have been developed into a sophisticated 
theory covering various aspects of dispersion forces, they can only render ap- 
proximate results valid in the non-retarded limit where the object separa- 
tions are sufficiently small so that the influence of the transverse electromag- 
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netic field can be disregarded. This was first demonstrated by Casimir and 
Polder [107,242]. Using a normal-mode expansion of the quantized electro- 
magnetic field inside a planar cavity bounded by perfectly conducting plates, 
they showed that the force between the plates can be derived from the total 
zero-point energy of the modes 

£ = E 5*"*- (6) 

k 

The difficulty that this energy is divergent was overcome by subtracting the re- 
spective diverging energy corresponding to infinite plate separation, the finite 
result implying a force per unit area 



7T 



He 1 



240 z 4 



(7) 



In a similar way, they obtained the force on an atom near one of such plates 
and the force between two atoms in free space from the ground-state energy 
of the respective system in leading-order perturbation theory. They recovered 
the results of the non-retarded limit, Eqs. (2) and (3), and found that in the 
retarded limit the atom-atom and atom-plate potentials are given by 

23fica(0)a / (0) 
U{Z) ~ 64^ (8) 

and 

_ 3fica(0) 

U[) 327r%o^ 4 ' [ ' 
respectively, which correspond to attractive 1/z 8 and 1/z 5 forces that decrease 
more rapidly than the ones in the non-retarded limit. Casimir and Polder had 
thus developed a unified theory to describe dispersion interactions over a large 
range of distances. 

Normal-mode techniques have since been widely used to study dispersion 
interactions. The two-atom interaction has been confirmed in various ways 
[243-257], inter alia by basing the calculations on the multipolar-coupling 
scheme [243,245,254] in place of the minimal-coupling scheme originally used 
by Casimir and Polder, and relativistic corrections have been considered [258]. 
Extensions include the interaction of three [247,255,259-262] or more atoms 
[263,264], the influence of higher-order multipole moments [265] and perma- 
nent dipole moments [266] on the two-atom force, the interaction between 
anisotropically polarizable atoms [267] and that between a polarizable and a 
magnetizable atom [268-274] (for an overview, cf. Ref. [275]). In particular, 
it was found that in the retarded limit the force between a polarizable and 
a magnetizable atom is repulsive as in the non-retarded limit, but follows 
the same 1/z 8 power law as that between two polarizable atoms. Further- 
more, the interaction of atoms in excited energy eigenstates [276-282] and 
the influence of external conditions such as finite temperature [273,283-286], 
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applied electromagnetic fields [285], or additional bodies [267,287-289] on the 
atom-atom interaction have been studied. In particular, when the interatomic 
separation exceeds the thermal wavelength, the force decreases more slowly 
(~ l/z 7 ) than in the zero-temperature limit. Similarly, the Casimir-Polder re- 
sult for the atom-plate interaction has been confirmed [247-250,252,269,290], 
atoms that carry permanent electric dipole moments [291] or are magnetizable 
[292] have been considered and the influence of finite temperature [284,293] 
as well as force fluctuations [294] has been studied. In close analogy to the 
atom-atom interaction, it was found that the interaction between a magne- 
tizable atom and a perfectly conducting plate is repulsive and that the force 
decreases more slowly (~ l/z 4 ) than in the zero-temperature limit as soon as 
the atom-plate separation exceeds the thermal wavelength. In contrast to the 
atom-atom interaction, the atom-plate potential for an atom in an excited 
energy eigenstate was found to show an oscillatory behavior in the retarded 
limit [251,287,295,296], thereby making the effect of the transverse electro- 
magnetic field more explicit. In addition, atoms interacting with bodies of 
different shapes and materials have been considered, such as: Perfectly con- 
ducting planar [291,297-301] and parabolic cavities [302,303]; metal [304], elec- 
tric [305-308] and magneto-electric half spaces [309]; electric planar [310,311] 
and spherical cavities [312]. 

Needless to say that the pioneering work of Casimir and Polder on disper- 
sion forces has also stimulated further studies of the problem of body-body 
interactions (for reviews see Refs. [44,313]). Apart from confirming and in- 
terpreting the original results [314-316], normal-mode techniques have been 
employed to include effects that arise from finite temperatures [317], surface 
roughness [318-320], the presence of a dielectric medium between the plates 
[321] and even virtual electron-positron pairs [322,323] (where the latter were 
found to be negligibly small). As in the case of atom-body interactions, vari- 
ous other geometries and materials have been considered such as: Two electric 
[324-327], dielectric [328], locally [329-331] and non-locally responding metal 
plates [332]; two plates that are polarizable and magnetizable [333-339]; the 
faces of a perfectly conducting rectangular cavity [340,341]; two electric multi- 
layer stacks [342]; a perfectly conducting plate and cylinder [343]; two electric 
spheres [344]; a perfectly conducting plate and a small electric sphere [345]; 
a sphere and a surrounding spherical cavity [346]. The results qualitatively 
resemble the findings for the atom-atom and atom-body interactions. In par- 
ticular, retardation was found to lead to a stronger asymptotic decrease of 
the forces which is softened due to thermal effects as soon as the separations 
exceed the thermal wavelengths; and the force between a polarizable object 
and a magnetizable one was found to be repulsive. Perhaps a more surpris- 
ing result is the fact that two birefringent plates may exert a non-vanishing 
dispersion torque on each other [347,348]. Moreover, the problem of Casimir 
energies of single bodies 12 has been addressed, motivated by a conjecture 

12 The Casimir energy of a single body can be defined as the geometry-dependent 
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made by Casimir [349], according to which an attractive Casimir energy of an 
electron (modeled as a small perfectly conducting sphere) should be able to 
counterbalance the repulsive self-energy of the electron charge and thus ex- 
plain its stability. 13 However, the energy of a perfectly conducting sphere was 
found to be repulsive [350,351], with similar findings for a weakly dielectric 
sphere [352-355]. On the contrary, the Casimir energy of a weakly dielectric 
cylinder was found to be attractive, in agreement with expectations [353,356]. 
The physical significance of Casimir energies of single objects is yet unclear; 
in particular it was shown by pairwise summation over microscopic dispersion 
energies that dispersion energies of macroscopic objects are in fact dominated 
by the always attractive volume and surface energies and may hence never be 
observed [357]. In standard calculations of Casimir energies, these volume and 
surface energies are either not considered from the very beginning or discarded 
during regularization procedures [355]. 

Normal-mode techniques have proved to be a powerful tool for studying disper- 
sion forces (cf. also Refs. [3,358]). Nevertheless, some principal limitations of 
the approach have become apparent recently, in particular in view of the new 
challenges in connection with recent improvements on the experimental side. 
So, normal-mode calculations can become extremely cumbersome when ap- 
plied to object geometries relevant to practice or when a realistic description of 
the electromagnetic properties of the interacting objects is required. The lim- 
itations are also illustrated by the controversy regarding the low-temperature 
behavior of dispersion forces on bodies (for a recent account of the debate, see 
Ref. [359] and references therein). The answer to this question requires detailed 
knowledge of the complicated interplay of positional, thermal and spectral fac- 
tors. To see this, one has to bear in mind that, in general, a large range of 
frequencies contributes to the forces where the relative influence of different 
frequency intervals is determined by the object-object separation, tempera- 
ture and the frequency dependence of the object properties. As a consequence, 
approximations such as long-/short-range, high-/low-temperature or perfect- 
reflectivity limits become intrinsically intertwined. A typical material property 
relevant to dispersion forces is the permittivity which is a complex function of 
frequency, with the real and the imaginary part being responsible for disper- 
sion and absorption, respectively. In particular, absorption which introduces 
additional noise into a system, inhibits the application of normal-mode expan- 
sion on a macroscopic level. This point was first taken into account by Lifshitz 
in his calculation of the dispersion force between two electric half spaces at 
finite temperature [360,361] where he derived the force from the average of the 
stress tensor of the fluctuating electromagnetic field at the surfaces of the half 
spaces, with the source of the field being the fluctuating noise current within 

part of the total electromagnetic energy where the notion geometry-dependent part 
is subject to ambiguities, cf. the discussion below. For further reading on the Casimir 
energy of a single body, cf. Ref. [313]. 
13 For a further a discussion of this idea, cf . Ref. [3] . 
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the dielectric matter. The required average was obtained by noting that the 
current fluctuations are linked to the imaginary part of the permittivity via 
the fluctuation-dissipation theorem. In this way, Lifshitz could express the 
force per unit area in terms of the permittivities s(u), s'(u) of the two half 
spaces where in particular, in the non-retarded (zero-temperature) limit the 
force per unit area, was obtained to be 



The Lifshitz theory has been applied and extended by a number of authors 
(for an overview, see Refs. [44,313]), who studied the influence of different 
frequency ranges [362-365], effects of finite temperatures [366-369] and sur- 
face roughness [370], and other planar structures such as electrolytic half 
spaces separated by a dielectric [371], magneto-electric half spaces [372], metal 
plates of finite thickness [373], metal half spaces exhibiting non-local proper- 
ties [368,374] and electric multilayer systems [375] (for further reading, cf., 
e.g., Ref. [233]). A typical approximation for treating small deviations from 
planar structures (like a sphere that is sufficiently close by a plate [330]) is 
the proximity force approximation where it is assumed that the interaction of 
two objects with gently curved surfaces can be obtained by simply integrating 
the (Lifshitz) force per unit area along the surfaces [376]. 14 While the debate 
regarding the temperature dependence of the force between realistic metal 
half spaces still seems unsettled, general consensus is reached that inclusion 
or neglect of material absorption (i.e., use of a Drude-type or a plasma- type 
permittivity) leads to the disagreeing results [359]. It is worth noting that 
the forces in a planar structure can be reexpressed in terms of (frequency- 
dependent) reflection coefficients directly accessible from experiments. This 
formulation of the theory has been applied to metal [378-380] and electric 
half spaces [381], metal half spaces with non-local properties [382], electric 
multilayer stacks [383,384] and, in some approximation, to rough perfectly 
conducting [385,386] and metal half spaces [387,388] where the surface rough- 
ness can give rise to a tangential force component [389] and a torque [390]. 

Lifshitz's idea of expressing dispersion forces in terms of response functions 
is of course not restricted to planar systems, but can be extended to arbi- 
trary geometries. This can be achieved by expressing the results obtained by 
normal-mode expansion in terms of the Green tensor of the (classical) elec- 
tromagnetic field [391-393]. Alternatively, the Green tensor which contains 
all the necessary information on the shape and the relevant electromagnetic 
properties of the objects, can be introduced by applying path-integral tech- 
niques [394,395] or employing the fluctuation-dissipation theorem [396]. The 



Recently, validity limits for the proximity force approximation in the case of per- 
fectly conducting objects have been discussed on the basis of numerical calculations 




(10) 



[377]. 
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theory has been used to study the forces between two perfectly conducting 
plates [391,397,398], a perfectly conducting plate and a perfectly permeable 
plate [393], two dielectric half spaces [328,399], two electric plates [396] and, in 
some approximation, two perfectly conducting spheres [393,397], a perfectly 
conducting sphere and a perfectly conducting plate [397,398]; the force on a 
rectangular piston [400]; and the force and the torque between weakly dielec- 
tric objects of arbitrary shapes [392]. In particular, it was shown that the 
force between two mirror-symmetric electric objects is always attractive [401]. 
Casimir energies of a perfectly conducting sphere [391], a dielectric sphere 
[395], a magneto-dielectric sphere [402] and a perfectly conducting cylinder 
[403] have also been studied in this way. 

The concept of linear-response theory has also been widely used to study 
dispersion forces on atoms. In particular, it can be shown that the (position- 
dependent part of the) interaction energy between a ground-state atom and the 
(body-assisted) electromagnetic vacuum in leading-order perturbation theory 
can be expressed in terms of the linear response functions of the two systems, 

= ^ rda 2 «(iOTrG«(r,r,iO, (11) 

i.e., the atomic polarizability a(u) on the one hand and the scattering Green 
tensor G^(r, r, u) of the electromagnetic field on the other, thus rendering 
a general expression for the force on an atom in the presence of arbitrary 
bodies [404-407] 15 (for an alternative, semiclassical approach based on find- 
ing the eigenenergies of the classical electromagnetic field interacting with a 
harmonic-oscillator atom, see Ref. [200]). The method which can easily be ex- 
tended to thermal fields [408-410], also applies to the dispersion interaction of 
two ground-state atoms in free space [404,408] or in the presence of bodies [216] 
(cf. also Refs. [200,411]). However, it cannot directly be applied to atoms in 
excited energy eigenstates where it is necessary to again start from the leading- 
order interaction energy and only express the field contribution in terms of the 
respective response function [412,413]. Linear response theory has been used 
to study the dispersion interaction of a single ground-state atom with a multi- 
tude of bodies such as: Perfectly conducting plates [404,405,407,408]; dielectric 
[414], electric [405,407,409,415] and magneto-electric half spaces [416]; metal 
half spaces exhibiting non-local properties [406,417] and/or surface roughness 
[406,418] or being covered by a thin overlayer [419]; perfectly conducting [420] 
and dielectric spheres [395,414,421,422]; dielectric cylinders [414,423-425]; per- 
fectly conducting and electric planar cavities [426]; dielectric spherical [427], 
cylindrical [424] and perfectly conducting wedge-shaped cavities [395,422]. 
Moreover, the force on an atom in an excited energy eigenstate in front of 

15 Note that the method is the natural extension of the approach based on the 
electrostatic response function which is now replaced by the response function for 
the complete electromagnetic field including its transverse part. 



16 



a perfectly conducting [412,413], dielectric [413] and birefringent dielectric 
half space [428] has been considered; and the interaction of two ground-state 
atoms embedded in a non-locally responding electrolyte [429], placed near a 
perfectly conducting [430] and a electric half space [415,431] or inside a per- 
fectly conducting [411,430] and dielectric planar cavity [432] has been studied. 

Finally, relations between microscopic and macroscopic dispersion forces have 
been established whose validity is no longer restricted to the non-retarded 
limit. Modeling macroscopic bodies as collections of harmonic-oscillator atoms 
interacting with the electromagnetic field and calculating the total energy 
of the interacting system, both the the force on a single ground-state atom 
in the presence of a dielectric half space and the force between two dielec- 
tric half spaces were derived from microscopic atom-atom interactions [433] 
where the former result was later extended beyond the harmonic-oscillator 
model [434]. A harmonic-oscillator model of atoms with the atoms being cou- 
pled to a heat bath, was used to derive the force between absorbing dielectric 
half spaces, confirming the result of Lifshitz theory [435]. The microscopic- 
model calculations show that only in the limit of weakly polarizable bodies, 
i.e., small values of the susceptibilities, a pairwise sum over two-atom forces 
is sufficient to obtain the total force, stressing once more the importance of 
many-atom interactions in the context of body-assisted dispersion forces. Pair- 
wise summation of two-atom interactions can be used to obtain an approx- 
imate description of the interaction between intricately shaped objects, e.g., 
bodies with rough surfaces [436], 16 or of atom-body/body-body interactions 
involving excited atoms and/or bodies [281,282,439]. Conversely, from well- 
known formulas for the body-body interaction, formulas for the atom-body 
interaction [360,361,394,395,440-446] as well as the atom-atom interaction 
[286,360,361,443,446,447] can be obtained in the limit of the respective sus- 
ceptibilities being asymptotically small. 

As we have seen, various concepts have been developed to describe disper- 
sion forces — an overview over the different scenarios which have been studied 
theoretically, is given in App. A in tabular form. These concepts, to some 
extent, are based on different basic assumptions and hence impose different 
limitations upon the applicability. The QED concepts based on normal-mode 
expansion of the quantized electromagnetic field typically suffer from the fact 
that when macroscopic bodies come into play, these bodies should be regarded 
as non-absorbing and hence also non-dispersing. To overcome this difficulty, 
arguments from other theories, such as the fluctuation-dissipation theorem 
of statistical physics, must be borrowed. On the contrary, the concepts based 
on linear response theory abandon an explicit field quantization and employ 



For bodies with small deviations from the planar geometry, the result of pairwise 
summation can be improved by introducing a correction factor obtained from Lifshitz 
theory [437,438]. 
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the fluctuation-dissipation theorem from the beginning. However, the appli- 
cability of methods that make use of the fluctuation-dissipation theorem by 
some means or other is limited to equilibrium systems — a disadvantage when 
dynamical aspects of excited atoms are to be considered. All concepts have in 
common that macroscopic bodies are typically described in terms of macro- 
scopic electrodynamics, i.e., boundary conditions at surfaces of discontinuity 
and/or constitutive relations. 

In this article we show that by following the formalism of macroscopic QED in 
media (as developed, e.g., in Refs. [448-450]) from the very beginning, one can 
obtain a unified approach to dispersion forces which does not only combine 
the benefits of normal-mode QED and linear-response theory in a natural way, 
but also accentuates the common origin of and intimate relations between the 
different types offerees and extends the range of application. In particular, the 
approach can be used to study dispersion forces for a wide class of different 
scenarios, including many of those listed in the above overview which have 
originally been studied by means of a variety of different methods. 

The further contents of the article are organized as follows. In Sec. 2, the 
main features of the quantization of the electromagnetic field in linear, dis- 
persing and absorbing media and the interaction of the medium-assisted field 
with atoms is outlined, with special emphasis on magneto-electric media de- 
scribed in terms of spatially varying permittivities and permeabilities which 
are complex functions of frequency. On this basis, in Sec. 3, very general for- 
mulas for the force on a macroscopic body due to its interaction with other 
macroscopic bodies are presented which are valid for arbitrarily shaped bod- 
ies as all the relevant properties of the bodies are fully expressed in terms 
of the Green tensor of the associated macroscopic Maxwell equations. Both 
Casimir stress and Casimir force are introduced, and a very general relation 
to many-atom van der Waals forces is established. In particular, it is shown 
that both the force on a single ground-state atom interacting with a body and 
the force between two ground-state atoms, can be obtained as limiting cases 
of the general formulas. In Sec. 4, forces on individual atoms interacting with 
the body-assisted electromagnetic field are studied in more detail, with special 
emphasis on explicitly solving the corresponding quantum-mechanical inter- 
action problem. It is demonstrated how the force on one or two ground-state 
atoms in the presence of magneto-electric bodies can be calculated by leading- 
order perturbation theory, the results agreeing with those obtained in Sec. 3. 
A number of examples is studied where it is shown that dispersion forces are 
often given by simple asymptotic power laws in the retarded and non-retarded 
limits. The force on a single atom initially prepared in an arbitrary excited 
quantum state is calculated by solving the atom-field dynamics, leading to 
explicitly time-dependent results. Some concluding are given in Sec. 5. 
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2 Elements of QED in linearly responding media 

It is well known that the properties of the electromagnetic field in media can 
significantly differ from those observed in free space, and hence, the interaction 
of the field with atoms can strongly be influenced by the presence of media. 
In classical electrodynamics, linear media are commonly described in terms of 
phenomenologically introduced macroscopic electric and magnetic susceptibil- 
ities (or permittivities and permeabilities, respectively) available from mea- 
surable data. This concept which can be transferred to quantum electrody- 
namics, has the benefit of being universally valid because it uses only very 
general physical properties, without the need for specific microscopic matter 
models and involved ab initio calculations. 

2.1 The medium- assisted electromagnetic field 

The medium-assisted electromagnetic field in the absence of free charges or 

currents obeys the macroscopic Maxwell equations which in the Fourier do- 
main read 

V-B(r,w)=0, (12) 

VxE(r,w) -icuB(r,cu) = 0, (13) 

e V-E(r,w)=p in (r,w), (14) 

K VxB(r,w) + iwe E(r,w) = j. n (r,cj) (15) 

{kq = fiQ 1 ) where the internal charge and current densities of the magneto- 
electric media p in (r, a;) and j in (r, a;) are the sources for the electric field E(r, a;) 

and the induction field B(r, uj). Note that the picture-independent Fourier 
components 0(t,u) of an operator field 0(r) are defined according to 

0(t)= duO(r,u) + H.c. (16) 
Jo 

so that 0(r,u,t) = e _lcJ ^ _t '^0(r, uj, f) in the Heisenberg picture. Since the 
internal charge and current densities are subject to the continuity equation 

-i^p. n (r,^) + V-j. n (r,^) = 0, (17) 

they may be related to polarization and magnetization fields P(r, a;) and 
M(r, a;) as follows: 

pJ r ,u)=-V-P(r,u), (18) 

j. (r, u)=- icuP(r, u) + V xM(r, u). (19) 
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Upon introducing the displacement field 



D(r, uj) = e E(r, id) + P(r, id) 



(20) 



and the magnetic field 



H(r, id) = /t B(r, id) - M(r, id), 



(21) 



the inhomogeneous Maxwell equations (14) and (15) can hence be written in 
the familiar equivalent form 



where the source terms associated with the internal charge and current den- 
sities are now contained in the displacement and magnetic fields. 

In particular in the case of linearly and locally responding magneto-electric 
media, Eqs. (20) and (21) take the form 



\k(t,u) =/i _1 (r,cj)] with e(r,u>) and /i(r, uj) being the (relative) electric per- 
mittivity and magnetic permeability of the media, respectively. Causality im- 
plies that e(r, uj) and fj,(r,u) which vary with space in general, are complex- 
valued functions of frequency with the Kramers-Kronig relations being satis- 
fied [451]. 17 According to the fluctuation-dissipation theorem, P N (r, cj) and 
M N (r, oj) are the (linear) noise polarization and magnetization, respectively, 
associated with the (linear) absorption described by the imaginary parts of 
e(r, a;) [Im e(r, a;) > 0] and /i(r, a;) [Im /x(r, a;) > 0]. For simplicity, in Eqs. (24) 
and (25) the material is assumed to be isotropic. 18 

Substituting Eqs. (20), (21), (24) and (25) into Eq. (23) and making use 
of Eq. (13), one can verify that the electric field obeys the inhomogeneous 



Note that both metals and dielectrics can be described in terms of their permit- 
tivity with the main difference being that the permittivity of a dielectric is analytic 
in the whole upper half of the complex frequency plane whereas that of a metal is 
commonly assumed to exhibit a simple pole at uj = [451]. 

18 The theory can be extended to arbitrary media, by starting from the general linear 
response relation between the current density and the electric field. Formulas in this 
article which do not explicitly refer to material properties (but solely via the Green 
tensor) are valid for arbitrary linear media [452]. 



V-D(r,w) = 0, 

V x H(r, id) + i^D(r, id) = 



(22) 
(23) 



P(r, id) = e„[e(r, id) - l]E(r, u) + P N (r, id), 
M(r, id) = «o[l - «(r, u)]B(r, id) + M N (r, id) 



(24) 
(25) 
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Helmholtz equation 



UJ 



VxK(r,u)Vx -e(r,u) 



E(r,w) = iujfx^ (r,uj) 



(26) 



where the source term is determined by the noise current density 

j N (r,w) = -iwP N (r lW ) + VxM N (r,u;). (27) 
Note that noise current density and noise charge density 

p N (r,u;) = -V-P N (r,u;) (28) 
fulfill the continuity equation 

-iujp N (r,uj) + V-\ N (r,uj) = (29) 
[recall Eqs. (17) (19)]. The solution to Eq. (26) can be given in the form 

E(r, uj) = kcvio / dV G(r, r', uj) -j^r', uj) (30) 

which, according to Eq. (13), implies that 

B(r,u)=iio J dVVxG(r,r' )W )-j N (r». (31) 

Here, G(r, r', a;) is the (classical) Green tensor which is defined by the equation 



UJ 



V x /c(r, uj) V x e(r, u;) 



G(r,r',u) = 6(t-t')I 



(I: unit tensor) together with the boundary condition 
G(r, r', a;) — > for |r — r'| — > 



oo. 



(32) 



(33) 



It should be pointed out that the Green tensor is uniquely defined by Eqs. (32) 
and (33) provided that the strict inequalities lme(r, uj) >0 and Im//(r, uj) >0 
hold. Note that it is an analytic function of uj in the upper complex half plane 
and and has the following useful properties (Ajj = Aji): 



G*(r,r» = G(r,r',-u;*) 
G(r,r',uj) = G T (r',r,w), 



(34) 
(35) 



J d 3 s |— Im/t(s,cj) 



V s xG(s,r,cu) • V s xG*(s,r'w) 



uj 2 1 
H — -Ime(s,a;) G(r, s, uj) ■ G*(s, r', uj)> = ImG(r, r', uj). (36) 
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Noise polarization and magnetization and hence noise current density, can be 
related to dynamical variables f A (r, uj) and f](r,cj) (A G {e, m}) of the system 
which consists of the electromagnetic field and the magneto-electric matter, in- 
cluding the dissipative system responsible for absorption, as follows [448,449]: 



P N (r,cu) = iy^Ime(r,w)f e (r,w), 



Hk, 



7Y 



h Im/i(r, uj) 



\| 7T/i |//(r, 



uj) 



(37) 



f m (r,cu) (38) 



with the f\(r, uj) and f|(r, u) being attributed to the collective Bosonic exci- 
tations of the system, 



fxi(r,u), /^(r'X) = 5 X \'5 u ,5(r - r')5(u - u'), 



f Xi (r,uj)Jx'i'(r',uj') 



= 0. 



(39) 
(40) 



By substituting Eqs. (37) and (38) into Eq. (30), on recalling Eq. (27), we may 
express the medium-assisted electric field in terms of the dynamical variables 
f\(r,ou) and f|(r, uj) to obtain 



E(r,w)= £ |dVG A (r,r>)-f A (r>) 



(41) 



X=e,m 



where 



u / h 

G e (r,r',w) = i—J Ime(r',w) G(r,r',u>), 



G m (r,r',u;) = i- 

c 



c<j / h 



ne 



ImK(r» [V'xG(r',r,cj)] T , 



(42) 
(43) 



so that, according to Eq. (16), 

E(r) = / dwE(r,cu) + H.c. 
Jo 



//"CO ^ 
dVy du;G A (r,r>).f A (r>) + H.c.. 

,\ < ./)/ " 



(44) 



Note that the relation (36) can be written in the more compact form 
I as u A (r,s,u;j-Li A (r,s,u;j = 

A=e,m * 



V /d^GA^s,^)^! 7 ^'^,^) = ^cu 2 ImG(r,r',cu). (45) 

J 71 



By starting from Eq. (41) and making use of the Maxwell equations in the 
Fourier domain, Eqs. (13) and (23) together with Eqs. (20), (21), (24), (25), 
(37) and (38), the other electromagnetic-field quantities such as B(r), D(r) 
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and H(r), can be expressed in terms of the dynamical variables f A (r, uj) and 
f|(r, uj) in a straightforward way. In particular, one derives 

B(r,u;) = - ]T / dW x G A (r, r', u) -f A (r', w), (46) 

A=e,m 

and hence, 



B(r) = /°°dcjB(r,cj) + H.c. 
Jo 

dV/ ^VxG A (r,r>)-f A (r',u;) + H.c. (47) 
Jo 1W 



A=e,m ' 



In view of the treatment of the interaction of the medium-assisted electro- 
magnetic field with atoms, it may be expedient to express the electric and 
induction fields in terms of potentials, 

E(r) = -Vy3(r)-A(r), (48) 
B(r) = VxA(r). (49) 

In Coulomb gauge, V- A(r) = 0, the first and second terms on the r.h.s. of 
Eq. (48) are equal to the longitudinal (||) and transverse (_L) parts of the 
electric field, respectively. From Eq. (44) it then follows that Vy3(r) and A(r) 
can be expressed in terms of the dynamical variables f A (r, uo) and fj[(r, uj) as 



r poo 

V0(r) = -E"(r) = - J2 /dVy dcJlG A (r,r>)-f A (r>) + H.c., (50) 
A(r)= £ /dV/ ^ ± G A (r,r',o;)-f A (r',a;) + H.c. (51) 

\=e,m J ° 1UJ 

Note that the longitudinal (transverse) part of a vector field is given by 

F HW(r) = / dV<5 ll(±) (r-r')-F(r') (52) 

where 

<5ll(r) = -w(-M, tf- L (r) = *(r)l-tfll(r) (53) 



,47rr 

(and for a tensor field accordingly). 

The commutation relations for the electromagnetic fields can be deduced from 
the commutation relations for the dynamical variables f A (r, u;) and f A (r, uj) as 
given by Eqs. (39) and (40). In particular, it can be shown [448,449] that the 
electric and induction fields obey the well-known (equal-time) commutation 
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relations 



^(r),4(r') 



-iHe^eu'kdkSfr - r'), 







A(r),^(rO 



(54) 
(55) 



Introducing the canonical momentum field associated with the (transverse) 
vector potential, 

ri(r) = -eoE-^r), (56) 
one can also prove that the (equal-time) commutation relations 



4(r),n i ,(r') 
A(r),iKO 



iM^(r-r'), 



ft* (r), 1^(0 



(57) 
(58) 



are fulfilled. 



It is an almost trivial consequence of the quantization scheme that upon choos- 
ing the Hamiltonian of the combined system to be 



#mf= ]T /d 3 r^ dw/jwfJ(r,w)-f A (r,w), 



A=e,m 

the Heisenberg equation of motion 

6 = 



h 



H mi , O 



generates the correct Maxwell equations in the time domain, 

VxE(r) + B(r) = 0, 
VxH(r) -D(r) = 0. 



(59) 



(60) 

(61) 
(62) 



Note that the Maxwell equations V-B(r) = and V-D(r) =0 are fulfilled by 
construction. 



The Hilbert space can be spanned by Fock states obtained in the usual way 
by repeated application of the creation operators f|(r, uj) on the ground state 
|{0}) which is defined by 

f A (r,cu)|{0}) = VA,r,cu. (63) 

Note that the ground state refers to the combined system of the electromag- 
netic field and the medium. In particular, from Eq. (41) together with the 
commutation relations (39) and (40) one derives, on using the integral rela- 
tion (45), 

({0}|4(r,cc04t(rV)|{0}> = Ti- l h^u 2 \mG lV {v,v\uj)5{uj-uj') ) (64) 
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which reveals that the ground-state fluctuations of the electric field are deter- 
mined by the imaginary part of the Green tensor — in full agreement with the 
fluctuation-dissipation theorem. 

It should be stressed that lme(r, u) > and Im/x(r, u) > are assumed to 
hold everywhere. Even in almost empty regions or regions where absorption 
is very small and can be neglected in practice, the imaginary parts of the 
permittivity and permeability must not be set equal to zero in the integrands 
of expressions of the type (44). To allow for empty-space regions, the limits 
Ime(r, uo) — > and Im /x(r, oo) ^0 may be performed a posteriori, i.e., after tak- 
ing the desired expectation values and having carried out all spatial integrals. 
In this sense, the theory provides the quantized electromagnetic field in the 
presence of an arbitrary arrangement of linear, causal magneto-electric bodies 
characterized by their permittivities and permeabilities where Ime(r,w) > 
and Im yu(r, uj) > 0. 

As outlined above, quantization of the electromagnetic field in the presence of 
dispersing and absorbing magneto-electric bodies can be performed by start- 
ing from the macroscopic Maxwell equations including noise terms associated 
with absorption, expressing the electromagnetic field in terms of these noise 
terms and relating them to Bosonic dynamical variables in an appropriate 
way (cf. also Refs. [453-455]). Alternatively, absorption can be accounted for 
by expressing the electromagnetic field in terms of auxiliary fields with the 
dynamics of the auxiliary fields being such that the Maxwell equations are 
fulfilled [456]. It has been shown that the two approaches are equivalent [457]. 
It is worth noting that the quantization scheme is in full agreement with the 
results of (quasi-) microscopic models of dielectric matter where the polariza- 
tion is modeled by harmonic-oscillator fields and damping is accounted for 
by introducing a bath of additional harmonic oscillators [458]. After a Fano 
diagonalization [459] of the total Hamiltonian of the system (which consists of 
the electromagnetic field, the polarization and the bath), an expression of the 
form (59) is obtained. A model of this type was first developed for homoge- 
neous dielectrics [458] and later extended to inhomogeneous dielectric bodies 
[460-462], including bodies exhibiting non-local properties [463]. In particular 
in the latter case, the differential equation (32) for the Green tensor obviously 
changes to an integro-differential equation and Eqs. (37) and (38) must be 
modified accordingly. 



2.2 Atom-field interaction 

Let us consider a system of nonrelativistic particles of masses m a and charges 
q a which form an atomic system, e.g., an atom or a molecule, interacting with 
the medium-assisted electromagnetic field. The Hamiltonian governing the 
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dynamics of the atomic system (briefly referred to as atom in the following) in 
the absence of the medium-assisted electromagnetic field is commonly given 
in the form 

#at = E P~ + \ I d 3 rp at (r)y3 at (r) (65) 

where p a t(r) and v?at(r), respectively, are the charge density and the scalar 
potential which are attributed to the atom, 

Pat(r) = J2^ 6 ( r ~ ( 66 ) 

a 

and the standard commutation relations 

[r a i,Pa'i'} = ift8aa>8ii>, (68) 
[r a i,r a 'i'} = = [p ai ,p a >i>} (69) 

hold. Obviously, (p&t(r) and p a t(r) obey the Poisson equation 

£ A^ at (r) = -p a t(r), (70) 

and the continuity equation 

iSat(r) + V-j a t(r) =0 (71) 
is fulfilled where the atomic current density j a t(r) reads 

Jat(r) = \ E Qa [rJ(r - r a ) + S(r - i a )i a ] . (72) 

a 

It may be useful [464,465] to introduce center-of-mass and relative coordinates 

r A = ^2 — f Q , r a = r a - r A (73) 

(rriA = J2a m a) with the associated momenta being 

Pa = E P«> Pa = - — Pa ( 74 ) 

a m A 

Combining Eqs. (65) and (74), the atomic Hamiltonian may be written in the 
form 

~ 2 —2 

H at = + y + \ I d 3 rp at (r)^ at (r) 
2m A ' 

~ + E^I«)H (75) 
Zm A 
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where E n and |n) are the eigenenergies and eigenstates of the internal Hamil- 
tonian. From the commutation relations (68) and (69) it then follows that the 
non- vanishing commutators of the new variables are 



[rAi,PAi'] = i^ii', 

r aii Pa'i 



■kx Ix m °' 
inOiii [ o na i 

m A 



(76) 
(77) 



In particular, when m a r /m A <C 1, then 



r aii Pa'i' 



(78) 



Further atomic quantities that will be of interest are the atomic polarization 
P at (r) and magnetization M at (r) [466], 



Pat(r) = y.q a r a [ da 5 (r - i A - ar a ) , (79) 
Jo v J 

rl r ■ ■ 

M«t(r) = |5Zg a / da a 8(r-r A -aT a ^ a xf a - F Q xT a 5(r-r A -<rt a ) 



(80) 

it can be shown that for neutral atoms, the atomic charge and current densities 
are related to the atomic polarization and magnetization according to 



Pat(r) = -V-Pat(r) 



j at (r) = P at (r) + VxM at (r) + j r (r) 



(81) 
(82) 

j r (r) = |Vx[p at (r)xr A -r A xP at (r)] (83) 

which is due to the center-of-mass motion, is known as the Rontgen current 
density [254,464]. Note that Eqs. (70) and (81) imply 



and 
where 



£ V^ at (r) = P at (r) 



(84) 



Expanding the delta functions in Eqs. (79) and (80) in powers of the rela- 
tive coordinates r Q , we see that the leading-order terms are the electric and 
magnetic dipole densities associated with the atom, 

P at (r) = d5(r-r A ), (85) 
M at (r) = iM(r - r A ) (86) 
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where the electric and magnetic atomic dipole moments read 

d = E QoFa = E 9ar Q , (87) 

a a 

m = |X]^aXr Q . (88) 

Note that the second equality in Eq. (87) only holds for neutral atoms. Using 
the atomic Hamiltonian (75) together with the commutation relation (77) and 
the definition (74), one can easily verify the useful relation 

E — (m\p a \n) = iu mn d mn (89) 
a m a 

\umn = {E rn — E n ) / k, d mn = (m\d\n)] which in turn implies the well-known 
sum rule 

1 g2 
77Z ^ ' ^mnjdnmdmn ~\~ d mn d nm ) ^ ' 7: /• (90) 

Zn m a zm a 



2.2.1 Minimal coupling 

Having established the Hamiltonians of the medium-assisted field and the 
atom, we next consider the atom-field interaction. According to the mini- 
mal coupling scheme (cf., e.g., Ref. [254]), this may be done by making the 
replacement p a h- > p a — q a A(r a ) in the atomic Hamiltonian (65), summing 
the Hamiltonians of the medium-assisted field and the atom and adding the 
Coulomb interaction of the atom with the medium-assisted field, leading to 
[448,465] 

H= E /d 3 r /* duhwfl(r,u)-f x (r,u) + j^m" 1 [p Q - q a A(r a ) 

X=e,m J J ° a 

+ \ y*d 3 rp at (r)^ at (r) + J d 3 r p at ( r M r ) 
= H mi + H at + H mt (91) 

where tp(r) and A(r) must be thought of as being expressed in terms of the 
dynamical variables f\(r, a>) and f|(r, u), according to Eqs. (50) and (51). 
Hence, the atom-field interaction energy reads 

#int = E Qaffia) ~ E ~ PaM*a) + E S~ A 2 (*«). (92) 
a a m a a ^ m a 

Note that the scalar product of p a and A(r a ) commutes in the Coulomb gauge 
used. 
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The total electromagnetic field in the presence of the atom reads 

5(r) = E(r)-V6*(r), S(r) = B(r), (93) 

X>(r) = D(r) - e Vy3 at (r), H(r) = H(r). (94) 

Obviously, B(r) and X>(r) obey the Maxwell equations 

V-B(r)=0, (95) 

V-X>(r)=p at (r), (96) 

and it is a straightforward calculation [448,465] to verify that the Hamiltonian 
(91) generates the remaining two Maxwell equations 

Vx£(r) +B(r) = 0, (97) 

VxW(r)-i>(r)=L(r) (98) 
and the Newton equations of motion for the charged particles, 

m a r a = q a £ (r a ) + \q a \r a x B(r Q ) - B(r a ) x r J (99) 



where 

p a -g a A(r a )l. (100) 



r a =m a 1 



In many cases of practical interest one may assume that the atom is small 
compared to the wavelength of the relevant electromagnetic field. It is hence 
useful to employ center-of-mass and relative coordinates [Eqs. (73) and (74)] 
and apply the long-wavelength approximation by performing a leading-order 
expansion of the interaction Hamiltonian (92) in terms of the relative particle 
coordinates r a . Considering a neutral atom and recalling Eq. (50), one finds 

H mt = -d.E"(f A ) PV Hi a) + E A- A 2 (?a). (101) 

a "<>a a £m a 

Note that the last term on the r.h.s. of Eq. (101) is independent of the relative 
particle coordinates and hence does not act on the internal state of the atom. 
When considering processes caused by strong resonant transitions between 
different internal states of the atom, it may therefore be neglected. 



2.2.2 Multipolar coupling 

An equivalent description of the atom-field interaction that is widely used is 
based on the multipolar-coupling Hamiltonian. 19 For a neutral atom, the tran- 



19 For an extension of the formulas given below to the case of more than one atoms, 
see Ref. [467]. 
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sition from the minimal-coupling Hamiltonian (91) to the multipolar-coupling 
Hamiltonian is a canonical transformation of the dynamical variables, corre- 
sponding to a unitary transformation with the transformation operator being 
given by 



U = exp 



d 3 rP at (r)-A(r) 



(102) 



where A(r) and P a t(r) are defined by Eqs. (51) and (79), respectively. This 
transformation is commonly known as the Power-Zienau-Woolley transfor- 
mation [468,469]; obviously it does not change f a , 

r' a = Urjfl = r a (103) 

and a straightforward calculation yields [448,465] 

V' a = Ut> a U ] = Pa - g a A(r a ) - J d 3 r S a (r) x B(r) (104) 

and 

f A (r,u;) = Uf x (T,u)fr = f A (r,c) + i-| dV P^r') • G A (r', r, u>) (105) 
where 

S«(r) = q a f a da a5 (r-r A - af c 



m, 
m A 



-y,qp*P [ daad(r - i A - crip) + — P at (r). (106) 

A a •'O 7 m A 



Now we may express the minimal-coupling Hamiltonian (91) in terms of the 
transformed variables to obtain the multipolar-coupling Hamiltonian in the 
form 



H= E |d 3 rrd^f;t(r lW ).f;(r )W ) + l/d 3 rP2(r) 

\=e, m J J ° Z£ J 

+ [Pa + J d 3 r3^(r)xB'(r)] 2 - / d 3 r P' at (r) -E'(r). (107) 

Here, E'(r) and B'(r), respectively, are given by Eqs. (44) and (47) with 
f\(r,iv) [f](r, u)] being replaced with f A (r, u) [f A (r, u)]. Note that r' a = r a , 
t' a = t a , t> A = T A , P at (r)=P at (r), 3' a (r) = 3 a (r), B'(r)=B(r), but 



E'(r) = E(r)+e 1 Pi- t (r) 



(108) 



which means that the transformed (medium-assisted) electric field E'(r) has 
the physical meaning of a displacement field, in contrast to E(r) which has 
the physical meaning of an electric field. 
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Hamiltonian (107) can be decomposed into three parts, 

H = H m{/ + H at > + H int > 



(109) 



where H m e is given by Eq. (59) with the primed variables in place of the 
unprimed ones, 



H w£ >= E j d 3 r J X 'duhwi;(r,u)-f' x (r,u), 

A=e,m 



(110) 



Hat' is the atomic Hamiltonian, 



f> /2 Ti' 2 i 



2mA 



12 



Pa 



+ S»'><»'i 



2mA 

and H in t' is the coupling term, 
H int , = -J d 3 rP' at (r)-E'(r) - J d 3 r M'Jv) -B'(r) 



(111) 



+ E^-[/d-r^(r)x6'(r) 



+ -/d 3 rp;.P' at (r)xB'(r) (112) 



where 



a &TTl a JO 

J113) 

Note that in contrast to the physical magnetization M at (r) [Eq. (80)], M at (r) 
is defined in terms of the canonically conjugated momenta rather than the 
velocities, as is required in a canonical formalism. The Hamiltonian (107) 
implies the relation 



m a r' a = p' a + d 3 rS' a (r)xB'(r) 



(114) 



and it is not difficult to see [recall Eqs. (100) and (104)] that m a r' a — m a r a . It 
should be pointed out that the eigenenergies E' n of the internal Hamiltonian 
in Eq. (Ill) may be different from the corresponding ones of the internal 
Hamiltonian in Eq. (75), because of the additional term contained in 



r) = l / dr/3 at (r)v4(r) + — / d 3 r 



2Sr 



(115) 



Accordingly, the eigenstates of the two internal Hamiltonians are not related 
to each other via the unitary transformation U [Eq. (102)] in general. 
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One of the advantages of the multipolar coupling scheme is the fact that it al- 
lows for a systematic expansion in terms of the electric and magnetic multipole 
moments of the atom. In particular, in the long- wavelength approximation, by 
retaining only the leading-order terms in the relative coordinates r' a , the in- 
teraction energy (112) reads 



H int ,= -d'-B'(r' A )-^.B'(r' A )+J2^- f' a xB'(v' A ) 



+ 



8mA 



d'xB'(r A ) +_ p ' A .d'xB'(r , A ) 



(116) 



where 



m 



q a 

2m r 



r'a * Pa- 



(117) 



Note that, in contrast to m [Eq. (88)], m is defined in terms of the canonical 
momenta. The first two terms on the r.h.s. of Eq. (116) represent electric and 
magnetic dipole interactions, respectively; the next two terms describe the 
(generalized) diamagnetic interaction; and the last term is the Rontgen inter- 
action due to the center-of-mass motion. For non-magnetic atoms, Eq. (117) 
reduces to the interaction Hamiltonian in electric-dipole approximation, 



H 



int' 



-d'-E'(f' 



+ Pk. d 'xB'(r' ) 

m A 



(118) 



which in cases where the influence of the center-of-mass motion on the atom- 
field interaction does not need to be taken into account, reduces to 



H 



int' 



-d'-E'(r A ). 



(119) 



3 Forces on bodies 



Electromagnetic forces are Lorentz forces. As known, the total Lorentz force 
F L acting on the matter contained in a volume V is given by 



d 3 r 



v 



p(r)E(r)+j(r)xB(r) 



(120) 



Here, the electromagnetic field acts on the the total charge and current den- 
sities p(r) and j(r), respectively, which in general include the internal charge 
and current densities p in (r) and jin(r), respectively, which are attributed to a 
medium [Eqs. (18) and (19)] as well as those due to the presence of additional 
sources, such as p at (r) and j a t(r) [Eqs. (66) and (72)]. With the help of the 
Maxwell equations [as given by Eqs. (12) (15) with pi n (r)i— > p(r), j in (r) i— >j(r)] 
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one easily finds 



p(r)E(r) + j(r) xB(r) = V- t(r) - e ^ 



E(r)xB(rl 



so that 



/ da- T(r) - e - , 

Jdv at Jv 



d 3 rE(r)xB(r) 



where the Maxwell stress tensor 

t(r) = 5 E(r)E(r) + ^- 1 B(r)B(r) - \ [e E 2 (r) + ^'(r) 



I 



(121) 
(122) 
(123) 



has been introduced. In particular, if the volume integral in the second term 
on the r.h.s. of Eq. (122) does not depend on time, then the total force reduces 
to the surface integral 



dF T 



dV 



(124) 



where 

dF L = da- t(r) = t(r)-da (125) 

may be regarded as the infinitesimal force element acting on an infinitesimal 
surface element da. Note that a constant term in the stress tensor does not 
contribute to the integral in Eq. (124) and can therefore be omitted. 



If the Minkowski stress tensor 

T (M) (r) = D(r)E(r) + H(r)B(r) - \ [£>(r)-E(r) + H(r)-B(r) 



/ 



= T(r) + P(r)E(r) - M(r)B(r) - \ [P(r)-E(r) - M(r)-B(r) / (126) 

(which agrees with Abraham's stress tensor [470]) is used in Eq. (124) [together 
with Eq. (125)] instead of the Maxwell stress tensor T(r) to calculate the force, 
one finds 

dF< M > = da-t( M V) = dF L 

+ da- (P(r)E(r) - M(r)B(r) - \ [p(r) -E(r) - M(r) -B(r)] /}, (127) 

and it is seen that in general 

dF L ^da-t (M) (r)- (128) 

That is to say, the use of the Minkowski stress tensor is expected not to 
yield the Lorentz force, in general. Indeed, a careful analysis and interpreta- 
tion of classical electromagnetic force experiments [471-478] shows that the 
(energy-momentum four-tensor associated with the) Lorentz force passes the 
theoretical and experimental tests and qualifies for a correct description of 
the energy-momentum properties of the electromagnetic field in macroscopic 
electrodynamics [479] (also see Sees. 3.1 and 3.2). 
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In classical electrodynamics, electrically neutral material bodies at zero tem- 
perature which do not carry a permanent polarization and/or magnetization 
are not subject to a Lorentz force in the absence of external electromagnetic 
fields. As already noted in Sec. 1.1, the situation changes in quantum elec- 
trodynamics, since the ground-state fluctuations of the body-assisted electro- 
magnetic field and the body's polarization/magnetization charge and current 
densities can give rise to a non-vanishing ground-state expectation value of 
the Lorentz force — the Casimir force [480] 



d 3 r {({0}| [p(r)E(r')+j(r)xB(r')] |{0}>} r ,_ (129) 

Here, the coincidence limit r' — > r must be performed in such a way that 
unphysical (divergent) self-force contributions are discarded after the vacuum 
expectation value has been calculated for r' ^ r, an explicit prescription will 
be given below Eq. (134). 

To calculate the Casimir force, let us consider linear media that locally respond 
to the electromagnetic field and can be characterized by a spatially varying 
complex permittivity e(r, ui) and a spatially varying complex permeability 
fj,(r,u). Following Sec. 2.1, we may write the medium-assisted electric and 
induction fields in the form of Eqs. (44) and (47). Provided that the volume of 
interest V does not contain any additional charges or currents, the charge and 
current densities that are subject to the Lorentz force (120) are the internal 
ones, p(r)=p in (r) and j(r) = j in (r) [recall Eqs. (18) and (19)]. Making use of 
Eqs. (24) and (25) together with Eqs. (28)-(32), one can easily see that 

p(r,u) = -s V-{[e(r,u) - l]E(r,w)} + M^V-j^r, u;) (130) 

and 

j(r,w)= - i(ve [e(r,u>) - l]E(r,w) 

+ Vx{« [l-«(r,o;)]B(r,a;)}+j N (r,a;). (131) 

Taking into account that E(r, a;) and B(r, a;) can be given in the forms (30) 
and (46), respectively, and that the Green tensor G(r, r',ui) obeys the differ- 
ential equation (32), one may perform Eqs. (130) and (131) to obtain 

p(r,u) = ^J dVV-G(r,r>)-j N (r',u;), (132) 



j(r,c)= jdV 



_ _ 

Vx Vx — - 



G(r,r>)-Ur',u;). (133) 



In this way, the fields p(r, u), j(r, u), E(r, u) [Eq. (30)] and B(r, u) [Eq. (31)] 
are expressed in terms of the noise current density j N (r, uS). Making use of 
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Eq. (27) together with Eqs. (37) and (38) and recalling the commutation rela- 
tions (39) and (40), one can easily calculate the ground-state correlation func- 
tion ({0}|j (r, u;)j N (r', u/)|{0}) which can then be used, on recalling Eqs. (30), 
(31), (132) and (133), to calculate all the correlation functions relevant to the 
Casimir force, as given by Eq. (129). The result is [481] 



r roc / fjj^ 

- / d 3 r / dcu — V-ImG(r, r',cu) 
ix Jv Jo \c 2 



+ Tr< / x 



VxVx — - 

c 1 



ImG(r, r', uS) x V' 



= --/ d 3 r f>(^V-G(r,r',iO 

ix Jv Jo \c z 



-TW /x 



Vx Vx + 



e 



G(r,r',iOxV 



(134) 



[(TrT)j = Tiji, V introduces differentiation to the left] where it is now appar- 
ent that in the coincidence limit r' — > r the Green tensor has to be replaced 
with its scattering part at each space point. In particular, when the mate- 
rial in the space region V is homogeneous, then the Green tensor therein can 
be globally decomposed into a bulk part G^(r, r', uj) and a scattering part 
G (1) (r,r», 



G(r,r» = G (0) (r,r» + G (1) (r,r', 



(reV). 



(135) 



In this case, the coincidence limit r' — > r simply means that the Green ten- 
sor G(r, r', uj) can be globally replaced by its well-behaved scattering part 
G (1) (r,r». 



According to Eqs. (123)— (125) , the Casimir force can be equivalently rewritten 
as a surface integral over a stress tensor [480], 



F = / da- T(r, rV_ 

JdV 



(136) 



where 



T(r,r') = ({0}|{e E(r)E(rO +/^ 1 B(r)B(r') 

-i[e o E(r).E(r')+^ 1 B(r).B(r')]/}|{0}> (137) 



which leads to 



T(r,r') = S(r,r')-|[TrS(r,r') 



(138) 
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with 



Sir 



du 



,0 = 5/ 

7T JO 

Pj roc 

= -- 

7T JO 



ImG(r, r', cj) - V x ImG(r, r', a;) x V' 



e 



G^r'.iO + VxG^r'.iOxV 



(139) 



Both Eq. (134) and Eq. (136) [together with Eqs. (138) and (139)] are valid 
for arbitrary bodies that linearly respond to the electromagnetic field, since 
the force is fully determined by the Green tensor of the classical, macroscopic 
Maxwell equations with the material properties entering the force formulas 
only via the Green tensor. Moreover, Eqs. (134) and (139) reveal that the 
force is proportional to h and hence represents a pure quantum effect. The 
results can be generalized to finite temperatures T in a straightforward way, 
by averaging in Eqs. (129) and (137) over the thermal state instead of the 
vacuum state. As a consequence, the r.h.s. of the first equalities of Eqs. (134) 
and (139) are modified according to [480] 



dto . . . 



dtd coth f — — — 
\2k B T / 



(140) 



(&b, Boltzmann constant) so that the final forms of these equations change as 



n poo 

- dCM) - 2AfeTX;(l-5<S n o)/(&») (141) 

71 J0 n=0 



with 



2%k B T 



(142) 



being the Matsubara frequencies. 

3. 1 Casimir stress in planar structures 



Let us apply the theory to a planar magneto-electric structure defined accord- 
ing to 



e(r, uo) = < 





u) 


z < 0, 


< e(u) 




< z < d, 




u) 


z > d, 


' H-(z 


u) 


z < 0, 


< fi{uj) 




< z < d, 






z > d 



(143) 



(144) 
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and restrict our attention to the zero-temperature limit. To determine the 
Casimir stress in the interspace < z < d, we need the scattering part of the 
Green tensor in Eq. (139) for both spatial arguments within the interspace 
(0 < z — z' < d). Since the component q of the wave vector parallel to the 
interfaces is conserved and the polarizations a = s,p decouple, the required 
Green tensor (as given in App. B) can be expressed in terms of reflection 
coefficients r a ± = r a ±(u>,q) (q — |q|) referring to reflection of waves at the 
right (+) and left (— ) wall, respectively, as seen from the interspace. Explicit 
(recurrence) expressions for the reflection coefficients are available if the walls 
are multi-slab magneto-electrics, cf. Eqs. (B.7) and (B.8). 20 In the simplest 
case of two homogeneous, semi-infinite half spaces, the coefficients r CT± reduce 
to the well-known Fresnel amplitudes, Eq. (B.10). 

In order to determine the Casimir force, it is clear for symmetry reasons that 
one requires the z component T zz (r) = T zz (r, r) r /_^ r of the stress tensor in the 
interspace < z < d which, upon using the Green tensor from App. B, can be 
given in the form [480] 

fy poo poo q 

T zz (r) = -— J o d^ dq±iJL(iOg(z,%,q) (145) 
where the function g(z,£,q) is defined by 



g{zM,q) = - 2 
- 2 



b 2 {l + n- 2 ) + q 2 {l 



n 



b 2 (l + n- 2 ) -q 2 (l -n" 2 ) 
+ (b 2 -q 2 )(l-n- 2 
-{b 2 -q 2 ){l-n~ 2 



p~ 2bd r r D~ l 
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e~ 2bd r^r^D' 1 



e~ 2bz r< 



_|_ e -2b(d-z) r 



D7 1 
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-2b(d-z), 



p+ 



d; 1 



(146) 



with 



n = n(iO = Ve(i^H^)> 



D a = D a (i£, q) = 1 - r a+ r a _e 



-2bd 



(147) 

(148) 
(149) 



According to Eq. (136), Eq. (145) [together with Eqs. (146)-(149)] gives the 
force per unit area between two arbitrary planar multilayer stacks of (locally 
responding) dispersing and absorbing magneto-electric material where the in- 
terspace between them may contain an additional magneto-electric medium. 
It is worth noting that many specific planar systems that can be addressed 
by means of Eq. (145) have been studied by using alternative methods: The 
most prominent example is the case of two electric half spaces separated by 
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For continuous wall profiles, Riccati-type equations have to be solved [482]. 
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vacuum, as first considered by Lifshitz [360,361] and later readdressed [362- 
365,378,381], inter alia based on normal-mode QED [325-327,329]. Extended 
scenarios range from electric half spaces separated by an electric medium (as 
studied by means of electrostatic theory [204,205,236,241], normal-mode QED 
[321] and an extended Lifshitz theory [371]) over electric plates of finite thick- 
ness (as addressed on the basis of the Lifshitz theory [373] and linear-response 
theory [396]), electric multilayer stacks (as treated by generalizing the results 
of the Lifshitz theory [375,383,384] as well as normal-mode QED [342]) to 
magneto-electric half spaces (as studied by means of Lifshitz theory [393] as 
well as normal-mode QED [333-337,372]). For purely electric multilayer sys- 
tems, various effects not being taken into account by Eq. (145), have also 
been addressed in a number of works, such as the influence on the force of fi- 
nite temperature [324,328,330-332,335,337,360,361,366-369,379,380,391], sur- 
face roughness [318-320,239,240,370,385-389] and non-locally responding ma- 
terials [208,332,368,374]. As outlined above, finite temperature can be easily 
included in Eq. (145) by applying Eqs. (141) and (142) [480], whereas surface 
roughness as well as non-local material response can be taken into account 
by returning to the more general formula (136) [together with Eqs. (138) and 
(139)] and specifying the Green tensor appropriately. 

To further (numerically) evaluate Eq. (145), knowledge of the £- and (/-depen- 
dence of the reflection coefficients which depend on the respective planar sys- 
tem (see, e.g., the examples studied in Refs. [445,446,483-485]), is required. 
Let us here restrict our attention to the limit of perfectly conducting surfaces, 
i.e., r p ± = —r s ± = 1 and assume that the wall separation d is sufficiently large, 
so that the permittivity and the permeability of the medium in the interspace 
can be replaced by their static values e= e(0) and \i= /i(0). It is then not 
difficult to calculate the simplified integrals in Eq. (145) analytically to obtain 
the attractive Casimir force per unit area, F— T zz (d) [recall Eq. (136)] which 
acts between two perfectly reflecting plates as 21 
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+ (150) 



240 V e \3 3efi J d 4 



In particular, when the interspace is empty (e — 1, /i — 1), then Eq. (150) 
reduces to 

F = W^' (151) 

in agreement with the famous result (7) first obtained by Casimir on the basis 



21 Note that the terms proportional to e~ 2b ( d ~ z ) ' m Eq. (146) give rise to divergent 
integrals in Eq. (145) in the limit z^> d which obviously results from the assumptions 
of frequency-independent response of the plates and the intervening medium and 
infinite lateral extension of the plates. Since in a realistic system these contributions 
are canceled by similar contributions occurring at the backs of the plates [480], they 
can be discarded. 



38 



of a normal-mode expansion [242] and subsequently rederived [315,317], inter 
alia based on classical orbits [393,397] or dimensional arguments [314]. 




In the same approximation, the force that acts on a perfectly conducting 
plate in a planar cavity bounded by perfectly conducting walls and filled with 
a magneto-electric medium obviously reads 

Wl -s?) (162) 

where d\ (d T ) is the left (right) plate-wall separation. On the contrary, use of 
the Minkowski stress tensor (126) leads, for fJ> — l, to [483] 

Comparing the two results for fx = 1, we see that |F| < |F( M )|. Introduction of 
a (polarizable) medium into the interspace between the plate and the cavity 
walls is obviously associated with some screening, thereby reducing the force 
acting on the plate. Since the internal charges and currents of the interspace 
medium are fully included only in the Lorentz force [recall Eqs. (120) (125)], 
the force based on the Minkowski stress tensor or an equivalent quantity un- 
derestimates the screening effect and is hence larger than the Lorentz force in 
general. 



3.2 Macro- and micro-objects 



Let us return to the general formula (134) and consider the Casimir force 
acting on electric matter of susceptibility x(r, u;) —e(r,u) — l in some partic- 
ular space region V in the presence of arbitrary linearly responding bodies 
(outside V) in more detail. If G(r, r', uj) and G(r, r', cu), respectively, denote 
the Green tensors in the absence and presence of the electric matter in V 
with both of them taking into account the bodies in the remaining space, the 
differential equation (32) for G(r, r',ui) can be converted into the Dyson-type 
integral equation 

G(r, r', u) = G(r, r', co) + ^J v d 3 s X (s, v)G(r, s, u) ■ G(s, r', u) (154) 

where, for r G V, the Green tensor G(r, r', uS) satisfies the same differential 
equation as the free-space Green tensor, 

G(r,r» = 5(r-r')/, (155) 



Vx Vx - 
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from which it follows that 

uj 1 

-V-G(r,r» = -V5(r-r') (156) 
c z 

for reV. Equations (154) and (156) imply (r G V, u real) 

V4mG(r,r» = -V-Im[x(r,u;)G(r, r',uj)}. (157) 
In a similar way, one finds that (r G V, uj real) 

-,2 



Vx Vx — - 



ImG(r, r', u)xV' = ^ Im[x(r, uj) G(r, r', w)] x V'. (158) 



Substituting Eqs. (157) and (158) into Eq. (134) one can then show that the 
Casimir force acting on an electric body of volume V\ which is an inner part 
of a larger electric body (occupying volume V) reads [481] 



F =-ir d 4{x, dv ^ i «» vTr i G<r - r ' i «]' 



-2/ da- X (r,iQ[G(r,r',iO] r , A (159) 



In particular, in the case of an isolated body, i.e., when the region V D V\ is 
empty apart from the electric matter contained in V\, then the surface integral 
can be dropped, hence 

F = - / d 3 r / dC- X (r, iO VTr[G(r, r', i^)]^. (160) 
In Jvi Jo c z 



Let us briefly compare Eq. (159) with the equation obtained on the basis of 
the Minkowski stress tensor, 

f( m) = t_f d 3 r r de e ^ ( i0] Tr[G(rj i0]r/ __ (161) 

In JVj Jo cr 

It differs from Eq. (159) by a surface integral, in general [481]. Hence, the two 
force formulas agree in the case of an isolated body where the surface integrals 
do not contribute to the force and both equations reduce to Eq. (160). In 
contrast to Eq. (159), application of Eq. (161) to any inner, homogeneous part 
of a body leads to the paradoxical result that the force identically vanishes, 
because of Vx(r, i£) = 0. In other words, the only atoms that are subject to a 
force are those at the surface of the body. On the contrary, it is known that 
the van der Waals forces on all atoms of a body contribute to the Casimir 
force (Sec. 3.2.1). 

We have seen that within the framework of macroscopic QED, Casimir forces 
on linearly responding bodies can be expressed in terms of the respective 
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Green tensor of the Maxwell equations for the body-assisted electromagnetic 
field. Hence, the main problem to be solved in practice is the determination 
of the Green tensors for the specific systems of interest. Since closed formulas 
for Green tensors are only available for highly symmetric systems (see, e.g., 
Ref. [482]), approximative and numerical methods are required. For exam- 
ple, one can start from an appropriately chosen Green tensor as zeroth-order 
approximation to the exact one and perform a Born expansion of the exact 
Green tensor by iteratively solving the corresponding Dyson-type equation. In 
particular, iteratively solving Eq. (154) yields the Born series 



K=l 



2K 



C J 



K 

J=1 JV 



x G(r, si,u;)-G(si,S2,u;) • • • G(s K ,r',u). (162) 



3.2.1 Weakly polarizable bodies, micro-objects and atoms 

The force formulas (159) and (160) which follow from macroscopic QED with- 
out involved microscopic considerations, do not only apply to electric macro- 
objects but also to micro-objects. Moreover they also allow for studying the 
limiting case of individual atoms and determining in this way even the disper- 
sion forces with which bodies act on atoms and atoms act on each other in the 
presence of bodies. To see this, let us consider dielectric bodies which may be 
typically thought of as consisting of distinguishable (electrically neutral but 
polarizable) micro-constituents (again briefly referred to as atoms), so that 
the Clausius-Mossotti relation [451,486] 

X(r,uj) = eo^^a^fl -J7(r)a(w)/(3eo)] _1 

= e -Vr)aH[l + xM)/3] (163) 

may be assumed to be valid where a{u) is the atomic polarizability and rj(r) 
their number density. 22 Let V\ be the volume of an isolated dielectric body 
of susceptibility x( r )^)- The Born series (162) and the Clausius-Mossotti 
relation (163) imply that when the body is sufficiently small and/or weakly 
polarizable, then the force, as given by Eq. (160), is essentially determined by 
the leading-order term proportional to x( r , ~£q ^(r^u;), so that Eq. (160) 
approximates to 



F = 



In JVi 



( d 3 r^(r) /°°da 2 «(iOVTrG«(r,r,iO (164) 
JVi Jo 



22 Note that Eq. (163) is consistent with the requirement that both a(u) and %(r, 
be Fourier transforms of response functions iff r/(r)a(0)/(3eo) < 1- 

23 For a small and/or weakly polarizable body that is an inner part of a larger body, 
see Refs. [443,481,487]. 
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where, according to the decomposition 

G(t,t',u) = G (0) (r,r',u;) + G (1) (r,r» (reV) (165) 

[cf. Eq. (135)], G^(r,r',u) is simply the scattering part of the Green tensor 
G(r, r',cu) of the system without the dielectric body under consideration. 

It can be easily seen that Eq. (164) may be rewritten as [481] 

F = f d 3 rr?(r)F(r) (166) 

where 

F(r) = -VC/(r) (167) 

with 

= r^ea(iOTrG^(r,r,iO (168) 

Z7T JO 

being nothing but the van der Waals potential of a single ground-state atom 
of polarizability a(u>) at position r in the presence of arbitrary linearly re- 
sponding bodies at zero temperature (Sec. 4). Note that in the limiting case 
when V\ — > and rj — > oo but rjVi = 1, such that V\ covers a single atom at 
position ta, then F reduces to F(r^) — the force acting on a single atom. Note 
that a relation of the kind (166) was already used by Lifshitz to deduce the 
dispersion force between a single atom and an electric half space from that 
between a dielectric and an electric half space [360,361]. 

Equation (166) reveals that the force acting on a weakly polarizable dielectric 
body is the sum of the forces acting on all body atoms due to their interaction 
with other bodies giving rise to the Green tensor G(r, r', uj). Let us consider 
in more detail the interaction of a weakly polarizable body with a second iso- 
lated dielectric body of volume V{ which is also weakly polarizable [x'(r,u) ~ 
Sq 1 r]'(r)a'(u)]. Denoting the Green tensor associated with all remaining bodies 
except for the two under consideration by G(r, r', a;), expanding G(r, r', u), by 
starting from G(r, r', uj) in the Born series [i.e., using Eq. (162) with G i— > G, 
Gh> G and V i— > V{\ and again omitting terms of higher than linear order in 
the susceptibility, Eq. (164) leads to [481] 

F= f d 3 r?7(r) / dV ?/(r)F(r, r') (169) 

where 

F(r,r') = -VC/(r,r'), (170) 

is the force with which an atom of polarizability a'iuj) at position r' acts on 
an atom of polarizability a(oo) at position r with 

C/(r,r') = -M rdU'a(iOa\iOTr[G(ry^O-G(r',r,iO] (171) 

Z7T JO L J 
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being the two-atom van der Waals potential [467,488]. Equation (169) clearly 
shows that the Casimir force is a volume force and not a surface force as 
could be suggested on the basis of the Minkowski stress tensor. According 
to Eq. (169), the Casimir force between weakly polarizable dielectric bodies 
is the sum of all two-atom van der Waals forces between the body atoms — 
a result which was already obtained by Lifshitz for the special case of two 
dielectric half spaces [360,361]. In fact, such a relation formed the basis of early 
calculations of dispersion forces between bodies [230,231] and it is still used for 
treating bodies exhibiting surface roughness [436] or containing excited media 
[281,282]. 

The force between a polarizable atom and a magnetizable one can be obtained 
in a similar way [489]. For this purpose, we consider the interaction of polar- 
izable atoms contained in the first, weakly polarizable body (volume V\) with 
a second, weakly magnetizable body of volume V{ and magnetic susceptibility 
£(r, uj) = /x(r, uj) — 1 = norf (r) (3' (u) where f3'(u) denotes the magnetizability of 
the atoms contained in V{. Again expanding the Green tensor G(r, r', uj) asso- 
ciated with all the bodies except for the first one by starting from the Green 
tensor G(r, r', uj) (associated with all the bodies except for the two under 
consideration) and retaining only the linear order in ((r,uj), one obtains 

G(r, r',uj) — G(r, r',uj) - d 3 s((s,uj) G(r,s,w)xV s ■V s xG(s,r',w). 

Jv{ 

(172) 

Upon substitution of Eq. (172), the force on the first body (164) can again 
be written as a sum over two-atom forces, Eqs. (169) and (170) where the 
potential of a polarizable atom interacting with a magnetizable one reads 
[489] 

2ir Jo 

x Tr{ [G(r, s, iO x V r ] • V s x G(s, r, i£) . (173) 



3.2.2 Many- atom van der Waals interactions 

In general, not only two-atom interactions but all many-atom interactions 
must be taken into account to obtain exact dispersion forces involving macro- 
scopic bodies. To illustrate this point, let us return to Eq. (168) and consider 
the interaction of a single atom with a dielectric body (volume V), whose sus- 
ceptibility x(r, w) is given by the Clausius-Mossotti relation (163). Recall from 
Sec. 3.2.1 that G(r, r', uj) denotes the Green tensor of the whole arrangement 
of bodies and G(r, r', uj) is the Green tensor of all (background) bodies except 
for the one under consideration. Substituting the Born series (162) [G ^ G, 
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G i— > G] into Eq. (168), one can write the atom-body potential in the form 

oo 

U(r) = £ Uk{v) (174) 

K=0 

where 



^o(r) 



2tt 



da 2 «0O^G w (r,r,iO 



(175) 



is the atomic potential due to the background bodies 24 and (K > 1) 

(-i) K hn 



U K (r) = 



2nc 2K 



1] id 3 sj X (sj,iO 
j=i 



'v 
x Tr 



/°°da 2X+2 «(iO 

JO 



G(r, Si, i£)- G(si, s 2 , iO • • • G(s^, r, i£) (176) 



is the contribution to the potential that is of Ki\\ order in the susceptibility 
of the dielectric body. To further treat the sum on the right-hand side of 
Eq. (176), the Green tensors therein are decomposed into singular and regular 
parts according to 



G(r,r',u) = 



](£)V-r')'+G'(r,r>) 



(177) 



and use is made of the Clausius-Mossotti relation (163). A somewhat lengthy 
calculation then leads to the result that [488,490] 



OO 1 

U(r) = U (r) +1:- 

K=l ^ ■ 



n fdhj^sj) 
j=i j 



U(r,s 1 , ...,s K ) 



(178) 



where 
C/(ri,r 2 , ...,r N ) 



(1 + S 2 n)^ JO 



/ da 2 V(iO---«iv(io 
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G'(ri, r 2 , iO • G'(r 2 , r 3 , i£) • • • G'(t n , n, if) (179) 



is the iV-atom van der Waals potential in the presence of arbitrary linearly 
responding (background) bodies at zero temperature and hence generalizes the 
free-space result given in Refs. [263,264]. In Eq. (179), the symbol S introduces 
the symmetrization prescription 



Sf(r 1 ,r 2 ,...,r N ) 



(2 - 5 2N )N U( _ p{N) 



£ /( r n(i),r n (2),---,r n (iv)) (180) 



where P(N) denotes the permutation group of the numbers 1, 2, . . . , N. Note 
that G'(r,r',u) = G(r,r',u) for r^r'. Clearly, for N = 2, Eq. (179) reduces 
to the two-atom potential already given, Eq. (171). 
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Note that Uq(t) = for free-space background, i.e., if there are no further bodies. 
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Equation (178) reveals that under very general conditions the interaction of a 
single ground-state atom with a dielectric body of Clausius-Mossotti suscepti- 
bility may be regarded as being due to all many-atom interactions of the atom 
in question with the body atoms. A relation of this type was first derived for 
the special case of a homogeneous half space filled with harmonic-oscillator 
atoms [174,175,433] and later extended to homogeneous dielectric bodies of 
arbitrary shapes with vacuum background by means of the Ewald-Oseen ex- 
tinction theorem [434]. Note that in close analogy to Eq. (178), the dispersion 
force between two dielectric bodies is due to all many-atom interactions of 
atoms in the first body with atoms in the second one; this was explicitly ver- 
ified for the cases of two homogeneous half spaces [175,433] and spheres [235] 
and was also shown for two homogeneous bodies of arbitrary shapes [234]. 



4 Forces on atoms 

As demonstrated in Sees. 3.2.1 and 3.2.2, forces on individual ground-state 
atoms in the presence of linearly responding bodies can be deduced from 
the forces between dielectric bodies of Clausius-Mossotti type in the limiting 
case of the bodies being weakly polarizable. Alternatively, these forces can be 
derived by explicitly studying the interaction of atoms with the body-assisted 
electromagnetic field according to Sec. 2.2. This approach to dispersion forces 
on atoms allows for studying the influence of the internal atomic dynamics on 
the forces; in particular, excited atoms can also be considered. 

4-1 Ground-state atoms 

Atoms initially prepared in their ground state will remain in this state provided 
that the body-assisted field is also initially prepared in the ground state. In 
addition, the atom-field interaction in this case will involve only virtual, off- 
resonant transitions of the atoms and the body-assisted field. Consequently, 
dispersion forces on ground-state atoms may adequately be described within 
the framework of time-independent perturbation theory. 

4-1.1 Single- atom force 

Following the idea of Casimir and Polder [107], one can derive the force on a 
single atom at zero temperature from the shift AE of the system's ground- 
state energy E for given center-of-mass position of the atom which arises from 
the atom-field coupling. The potential U(ta), whose negative gradient gives 
the sought van der Waals force, is the position-dependent part of this energy 
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shift, 25 

AE = A {0) E + U(r A ). (181) 
It can be seen that the effective center-of-mass Hamiltonian 

H eS = ^ + U(r A ) (182) 
lm A 



leads to the equation of motion 

] 



F(r A ) = m A r A = -— [# eff , [# eff , m A r A ]] = -V A U(r A ). (183) 



Making use of the interaction Hamiltonian (119), one obtains the ground- 
state energy shift for a non-magnetic atom in second-order (i.e., leading-order) 
perturbation theory 26 as 

Ap vv [a* -dC" a |(0|({0}|-d-E(rA)|lA(r,^))|A;)| 2 

AE = > > d^rP du — r (184) 

[P, principal part; |l A (r, u)) = fj(r, w)|{0})]. Using Eqs. (39)-(41), (44) and 
(45), one derives 

AE = ^J2V f°° -^—u 2 d 0k -lmG(r A ,r A ,u;)-d k0 (185) 

7T ^ JO U k o + UJ 

where G(r, r', uj) is the Green tensor of the body configuration considered. 
By discarding the position-independent contribution AE^ associated with 
G^°\r A ,r A ,u) [recall Eqs. (135) and (181)] which may be thought of as being 
already included in the unperturbed energy, the van der Waals potential can 
be written in the form [465,491,492] 

U(r A ) = ^° | o °°da 2 Tr[a(iO-G«(r A ,r A ,iO] (186) 

where 

Q(w) = Sft?^^7^ (187) 



25 The position-independent part A^E is a contribution to the Lamb shift in free 
space; for a discussion of the Lamb shift within the multipolar coupling scheme, 
see, e.g., Ref. [468]. 

26 In the following, the multipolar coupling scheme will be employed and the primes 
discriminating the respective atomic and field variables from the minimal coupling 
ones will be dropped, for notational convenience. Equation (119) can be employed, 
because the second term in Eq. (118) gives rise to a contribution of order v/c (v, 
center-of-mass speed) [465] which can be neglected. 
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is the ground-state polarizability of the atom. For isotropic atoms, it simplifies 
to 

a(w) = a{u)l = hm A E / fco|d f |2 /, (188) 
so the potential simplifies to 

U(r A ) = ^ rda^OOTrG^^,^,^). (189) 
2n Jo 

The perturbative result hence agrees with what has been inferred from the 
force on weakly polarizable bodies and renders an explicit expression for the 
polarizability. Note that the scattering Green tensor G^ in Eq. (189) has 
exactly the same meaning as G in Eq. (168). From Eq. (189) together with 
Eq. (188) it can be seen that the potential can be given in the equivalent form 

U(r A ) = -^ I" duu 2 lm\a{uj)TrG {1) {r A ,r A ,uj)} (190) 
2tt Jo 1 J 

which allows for a simple physical interpretation of the force as being due to 
correlation of the fluctuating electromagnetic field with the corresponding in- 
duced electric dipole of the atom plus the correlation of the fluctuating electric 
dipole with its induced electric field [409]. It should be pointed out that an 
analogous treatment based on the minimal-coupling Hamiltonian (101) leads 
to the formally same result [465] where of course the unperturbed eigenstates 
and energies occurring in the polarizability (187) are now determined by the 
atomic Hamiltonian (65) in place of (111). Needless to say that both results 
are approximations to the same Hamiltonian of the total system. Bearing in 
mind that the ground-state energy shift is entirely due to virtual, off-resonant 
transitions, it is crucial to retain the A 2 term in the minimal coupling scheme 
which contributes to the ground-state energy shift already in first-order per- 
turbation theory. 

Equation (186) gives the potential of a single ground-state atom in the presence 
of an arbitrary arrangement of linearly responding bodies at zero temperature 
in terms of the polarizability of the atom and the Green tensor of the body- 
assisted electromagnetic field. A relation of this kind was first derived for arbi- 
trary electric bodies on the basis of linear-response theory [404,405,407], recall 
Eq. (11); alternatively, it was obtained from a QED path-integral approach 
[394] and semiclassical considerations [200]. A perturbative derivation based 
on macroscopic QED very similar to that presented here is given in Refs. [493- 
495]. The linear-response approach has also been applied to magneto-electric 
bodies [416] and finite temperatures [408-410]. Note that in close analogy to 
the case of the Casimir force, the single-atom potential at finite temperature 
can be obtained from the zero-temperature result (186) or (190) by making 
the replacements (141) or (140), respectively. 
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4-1.2 Two-atom force 



The interaction potential of two polarizable ground-state atoms in the pres- 
ence of linearly responding bodies giving rise to the Green tensor G(r, r',cu) 
can also be obtained by means of time-independent perturbation theory. The 
leading contribution is now of fourth order in the atom-field interaction and 
a somewhat lengthy calculation yields [467,488] 

x T±[a A (iO-G(r A ,r B ,%)-a B (%)-G(r B ,r A ,%)] (191) 
which for isotropic atoms reduces to 

U(r A ,r B ) = /"daVOO^OO 

x Tr[G(r A ,T B ,%)-G(r B ,r A ,%)]. (192) 

Note that Eq. (192) agrees with Eq. (171) [G'(ri,r 2 ,w) i-> G(r A ,r B ,cj) for 
ri^r 2 ]. The total force acting on atom A (B) can be obtained by supple- 
menting the single- atom force F(r A (B)) with the two- atom force 

F(T A {B),r B (A)) = -V A (B)U (r A ,r B ) (193) 

where in general F(r A ,r B ) ^ — F(r B ,r A ), due to the presence of the bodies. 

Equations (191) (193) form a general basis for calculating two-atom poten- 
tials in the presence of linearly responding bodies at zero temperature. An 
equation of the type (192) was first derived for the case of electric bodies 
by means of linear-response theory [216]. Derivations based on semiclassical 
models of harmonic-oscillator atoms interacting in the presence of perfectly 
conducting [411] and electric bodies [200] were given and extended to allow 
for finite temperatures [430]. Equation (192) has been used to study the inter- 
action of two atoms embedded in an electrolyte [429], situated near perfectly 
conducting [430], non-local metallic [215], electric [216,415,431] and magneto- 
electric half spaces [467,489], placed inside perfectly conducting [411,430] and 
electric planar cavities [432]. 

Let us consider the simplest case of two atoms in free space where the single- 
atom force identically vanishes and the two-atom potential (192) can be calcu- 
lated by using the free-space Green tensor G(r, r', oj) = Gf rce (r, r', a;) (App. B), 
leading to (r = (r^ — r B \) 

u ( r A,T B ) = - 3 dtaA(iZ)a B mg(tr/c) (194) 
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where 

g(x) = 2e~ 2x (3 + 6x + 5x 2 + 2x 3 + x 4 ) , (195) 

which shows that the interaction between two polarizable ground-state atoms 
is always attractive. Equations (194) and (195) are in agreement with the 
famous result of Casimir and Polder [107], whose derivation was based on 
fourth-order perturbation theory and normal-mode QED. The problem has 
been reconsidered many times, inter alia within the frameworks of normal- 
mode QED [243,244,246,250-252,256,257] and linear-response theory [404]; 
and it has even become a common textbook example [245,253,254]. 

In the retarded limit where r ^> c/u min (u m i n denoting the minimum of the 
relevant resonance frequencies of atoms A and B), the function g(£r /c) effec- 
tively limits the £ integral in Eq. (194) to a range where «A(B)(iO — Q; a(b)(0), 
so the potential approaches 

U(r A ,r B ) = M^p—> (196) 

as already pointed out in Ref. [107], cf. Eq. (8). In the non-retarded limit where 
r <C c/uj max (w max denoting the maximum of the relevant resonance frequencies 
of atoms A and B), the integral is limited by the polarizabilities cnA(B){i£) to 
a range where <?(i£) — g(0) = 6, leading to 

3/1 f°° 

U{Ta ^ b) = "16^6 I dWi£We (197) 

Upon recalling Eq. (188), one may easily verify that this non-retarded asymp- 
tote is nothing but the well-known London potential (2) which was origi- 
nally obtained from a perturbative treatment of the Coulomb interaction [164] 
(cf. Refs. [169,170] for similar derivations). 

It is illustrative to compare the potential between two polarizable atoms with 
the potential between a polarizable atom A of polarizability a A (u) and a 
magnetizable atom B of magnetizability Pb(u) which, according to Eq. (173), 
is given by 

U(r A ,r B ) = -M rd^ 2 »A OOMiO 

Z7T Jo 

xTr{[G(r A ,r,iOxVl-VxG(r,r A ,iO} . (198) 

1. L J J r=r B 

When the two atoms are in free space so that G(r, r', a;) = Gf rcc (r, r', a;), then 
Eq. (198) reads 

U(r A , r B ) = jT d £ e^ A mMmir/c) (199) 
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where 

h(x) = 2e~ 2x (l + 2x + x 2 ) (200) 

which is in agreement with results found on the basis of normal-mode QED 
[270,274]. In contrast to the attractive interaction between two polarizable 
atoms, the interaction between a polarizable and a magnetizable atom is al- 
ways repulsive, as can be easily seen from Eq. (199) together with Eq. (200). 
In particular in the retarded and non-retarded limits, respectively, Eq. (199) 
reduces to 

v(iAiB) = ^Mm (201) 

and 

U(r A ,r B ) = ^ f d ^ 2 a A (iO(3 B (iO (202) 
which was already given in Refs. [268,269] and [178,180]. 

Comparing Eqs. (196) and (201), we see that in the retarded limit the ab- 
solute value of the force between two polarizable atoms and that between a 
polarizable and a magnetizable atom follow the same 1/r 8 power law with 
the strength being weaker in the latter case by a factor 7/23. Comparison of 
Eqs. (197) and (202) shows that in the non-retarded limit the absolute value 
of the force between a polarizable and a magnetizable atom which follows a 
1/r 5 power law, is more weakly diverging than that between two polarizable 
atoms which follows a 1/r 7 power law. 



4-1.3 Atom in a planar structure 

Let us return to Eq. (186) for the potential of a single ground-state atom 
in the presence of an arbitrary arrangement of linearly responding bodies 
at zero temperature. It has been used to calculate the potential for a va- 
riety of particular geometries, including different planar structures (see be- 
low) as well as perfectly conducting [221,420], non-local metallic [222,223], 
dielectric [221,395,414,421,422,496] and magneto-electric spheres [488]; dielec- 
tric [226,414,423,424] and non-local metallic cylinders [224,225]; magneto- 
dielectric rings [490]; electric cylindrical shells [493-495]; electric [427] and 
non-local metallic spherical cavities [229]; electric cylindrical cavities [424] 
and perfectly conducting wedge-shaped cavities [395,422]. 

To illustrate the theory, let us consider an atom placed in a free-space region 
between two planar magneto-electric walls, as defined by Eqs. (143) and (144) 
with e(uj) = ji{uj) = 1. Inserting the Green tensor for this system (App. B) into 
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Eq. (189) leads to the single-atom potential [491,497] 



U(z A ) 



8tt 2 



dq r 



-2bz A 



+ e 



-2b(d-z A ) 



q 2 c 2 ^ 



e 



<_p+ 



(203) 



where z A is the separation of the atom from the left wall, d is the separation of 
the two walls, b and D a are given by Eqs. (148) and (149), respectively, with 
n(i£) = 1, and r a ± = r a ±(£,q) are again the reflection coefficients associated 
with the left/right walls. If the atom is placed near a single wall, say the right 
wall is missing, Eq. (203) reduces to (r a = r a -) 



-2bz A 



r s - 1 + 2 



q 2 c 2 " 



(204) 



4.1.3.1 Perfectly reflecting plate Consider first an atom placed near 
a perfectly reflecting electric (i.e., perfectly conducting) plate, r p = —r s = 1. 
By changing the integration variable in Eq. (204) from q to b [Eq. (148)], the 
resulting integral can be performed to obtain 



U{z A ) 



h 



1Qtv 2 £ z a 



/ d£a(ifle- 
Jo 



2iz A /c 



1 + 2 



+ 2 



, (205) 



in agreement with the result found by Casimir and Polder on the basis of 
normal-mode QED [107] (cf. also Refs. [247,250]) which has been reproduced 
by means of linear-response theory [404,405,407] and dynamical image-dipole 
treatments [290]. In the retarded limit, z A ^> c/u min , the exponential e~ 2 ^ ZA ^ c 
effectively limits the £ integral to the region where cc(i£) — «(0), so the poten- 
tial approaches 

U{z A )= 3hca{0) 



327r 2 e z A 



(206) 



as already noted in Refs. [107,245,249], cf. Eq. (9). In the non-retarded limit, 
z A <^ c/cj max , the polarizability a(i£) restricts the integration to the region 
where £z A /c~0, leading to 



U(z A ) = 



48tt£o4 k 



l d ofc| 2 = 



|d 2 |0) 



487T£o^A 



3 ' 



(207) 



in agreement with the result (3) obtained by Lennard- Jones on the basis of 
an electrostatic calculation [177] [recall Eq. (188)]. 



In contrast, the potential of an atom in front of an infinitely permeable plate 
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Fig. 1. The image dipole construction for an (a) electric (b) magnetic dipole in front 
of a perfectly reflecting electric plate, is shown. 



is repulsive, as can be seen by setting r p = —r s = — 1 in Eq. (204), 
h 



U(z A ) 



167r 2 e z\ Jo 
It approaches 



df a(i£)e 



-2£z A /c 



1 + 2 ^)+2 ^ 



U{z A ) 



3hca(0) 
32n 2 e z\ 

in the retarded limit (cf. also Ref. [309]) and 



U(z A ) 



(0|d 2 |0) 

48lTEoZ A 



(208) 



(209) 



(210) 



in the non-retarded limit. 



The different signs of the non-retarded potentials (207) and (210) in the cases 
of a perfectly reflecting electric and a perfectly reflecting magnetic plate, re- 
spectively, can be understood from an image-dipole model [177]. The non- 
retarded potential can be regarded as being due to the Coulomb interaction of 
an electric dipole d= (d x ,d y ,d z ) situated at distance z A from the plate with 
its image d' = (—d x ,—d y ,d z ) in the plate [Fig 1(a)]. The average interaction 
energy of the dipole and its image hence reads [451] 27 



U(z A ) = 



1 (0|d-d'-344|0) 



\d 2 + d 2 z \0} (0|d 2 |0) 



(211) 



2 4:7!-eo(2z A ) 3 647r£ -2i 48ne z 3 A 

[(0\d 2 x \0) = (0\d 2 y \0) = (0\d 2 z \0) = (l/3)(0|d 2 |0)], in agreement with Eq. (207). 

Since the interaction of a polarizable atom with an infinitely permeable plate is 
equivalent to the interaction of a magnetizable atom with a perfectly reflecting 
electric plate by virtue of the duality of electric and magnetic fields, a magnetic 



27 The factor 1/2 in Eq. (211) accounts for the fact that the second dipole is induced 
by the first one. 
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dipole rh= (m x ,m y ,rh z ) in front of a perfectly reflecting electric plate can be 
considered. A magnetic dipole behaves like a pseudo-vector under reflection, 
so its image is given by m'= (rh x ,rh y , —rh z ) [Fig 1(b)]. The interaction energy 
of the magnetic dipole and its image reads 

1 (0|m-m'-3m z m;|0) (0|m 2 |0) 
U ^)- 2 Ane (2z A r ~48ne z\ 

which by means of a duality transformation is equivalent to Eq. (210). The 
different signs of the potentials (207) and (210) can thus be understood from 
the different reflection behavior of electric and magnetic dipoles. 



4.1.3.2 Magneto-electric half space To be more realistic, let us next 
consider an atom in front of a semi-infinite half space of given s{uj) and fJ>(co). 
Upon substitution of the Fresnel reflection coefficients (B.10), Eq. (204) takes 
the form [491,497] 



b 



At(iO& + &i 
q 2 c 2 " 



-1 + 2 



(213) 



with bi = b l _ defined as in Eq. (B.9), in agreement with the result found 
by means of linear-response theory [416]. From Eq. (213), the results ob- 
tained by means of normal-mode QED [252,306,309] and linear-response the- 
ory [216,405,407] for an electric half space can also be recovered. 

One can show that in the retarded limit z A c/cj min (with u m i n being the 
minimum of all relevant atom and medium resonance frequencies) the potential 
takes the asymptotic form [491,497] 28 



U(z A ) = 



3hca(0) 

64lT 2 E Z A 



s: 



dv 



'(- 




\v 2 


v 4 ) 



e(0)v - yje(0)fi(0) - 1 + v 2 



1 ti(0)v - ^(0)/i(0) 



(214) 



which can be attractive or repulsive, depending on the strengths of the com- 
peting magnetic and electric properties of the half space. Figure 2 shows the 
borderline between attractive and repulsive potentials in the e(0)/i(0)-plane. 



28 Obviously this limit does not apply for metals where the condition can never be 
fulfilled. Similarly, Eqs. (215), (220) and (221) only hold for dielectrics. 
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Fig. 2. Borderline of attractive and repulsive retarded potentials of a ground-state 
atom in front of a magneto-dielectric half space. 

In particular, repulsion occurs iff /x(0) — 1 > 3.29[e(0) — 1] or /i(0) > 5.11e(0) for 
weak and strong magneto-dielectric properties, respectively. 

In the non-retarded limit, ti(0)za <C c/oj max [oj max , maximum of all relevant 



atom and medium resonance frequencies; n(0) = ye(0)/i(0)], the situation is 
more complex, because electric and magnetic medium properties give rise to 
potentials with different asymptotic power laws. In particular, the potential 
approaches 



in the case of a purely dielectric half space — in agreement with the result (4) 
found by considering the Coulomb interaction and using image-dipole methods 
[194] or linear-response theory [199] — and 



in the case of a purely magnetic one [491,497]. Thus for a magneto-dielectric 
half space the attractive l/z\ potential associated with the polarizability of 
the half space will always dominate the repulsive 1/za potential related to its 
magnetizability. It should be noted that the non-retarded limit is in general 
incompatible with the limit of perfect reflectivity considered in Sec. 4.1.3.1. 
So, Eq. (216) does not approach Eq. (210) as tends to infinity. On the 
contrary, Eq. (215) converges to Eq. (207) as tends to infinity [304]. 

Combining the observations from the retarded and non-retarded limits, one 
may thus expect a potential barrier, provided that the permeability of the 





(215) 




(216) 
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Fig. 3. The potential of a ground-state two-level atom in front of a magneto-dielectric 
half space, Eq. (213), is shown as a function of the distance between the atom and the 
half space for different values of /x(0) (ujp e /uJiQ = 0.75, wxe/^io = 1-03, tuxm/wio = 1, 

7eMo = 7mMo = 0.001). 

half space is sufficiently large [491,497]. This is illustrated in Fig. 3 where 
the potential of a two-level atom as a function of its distance from the half 
space, is shown for various values of the (static) permeability. In the figure, 
the permittivity and permeability of the half space have been assumed to be 
of Drude-Lorentz type, 



e(u) = 1 + 



n(ou) = 1 + 



UJ 



Pm 



(217) 



From the figure it is seen that with increasing value of /x(0), a potential barrier 
begins to form at intermediate distances, as expected. It is shifted to smaller 
distances and increases in height as the value of /x(0) is further increased. 



4.1.3.3 Plate of finite thickness Consider now an atom in front of 
magneto-electric plate of finite thickness d. Evaluating the relevant reflection 
coefficients (B.7) and (B.8) (d = d}_), one finds that Eq. (204) takes the form 
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duj w /c = 0.01 
dw w /c = 0.10 
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Fig. 4. The potential of a ground-state two-level atom in front of a magneto-dielectric 
plate, Eq. (218), is shown as a function of the atom-plate separation for different 
values of the plate thickness d [/x(0) = 5; whereas all other parameters are the same 
as in Fig. 2]. 



[491,497] 

[fi 2 {iOb 2 -b\] tanh(M) 



x 



2//(if )6&i + [fi 2 (iOb 2 + b\] tanh(M) 

[£ 2 (iOfe 2 -&?] tanh(M) 




2e(iO&6i + [e 2 (iO& 2 + 6?] tanh(M) 



(218) 



which reduces to the result in Ref. [310] in the special case of an electric plate. 
For an asymptotically thick plate, d^$> za, the exponential factor restricts the 
integral in Eq. (218) to a region where bid ~ dj (2za) 1. One may hence make 
the approximation tanh(&id)~l, leading back to Eq. (213) which demonstrates 
that the semi-infinite half space is a good model provided that d^$> za- 

On the contrary, in the limit of an asymptotically thin plate, fi(0)dCZi, the 
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approximation frid-Cl results in [491,497] 
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which in the retarded limit approaches 



U(z A ) 



hca(0)d 
160tt 2 Soz a 



14e 2 (0)-9 6^ 2 (0)-l 



s(0) 



(220) 



In the non-retarded limit one can again distinguish between a purely dielectric 
plate and a purely magnetic plate. Equation (219) approaches 



U(z A ) = 
in the former case and 

TTl \ M 

U{z A ) = —r- 2 2 , 

64ir 2 e z A Jo 



3hd 



64ir 2 e z A 
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(221) 



^ (-) «(iO 



[/i(iQ-l][3MiO + l] 



(222) 



in the latter case. Comparing Eqs. (220)-(222) with Eqs. (214)-(216), we see 
that the power laws change from zj^ to z A n+1 ^ when the plate thickness 
changes from being infinitely large to being infinitely small. 

If the permeability is sufficiently big, a magneto-dielectric plate of finite thick- 
ness features a potential wall [491,497], as illustrated in Fig. 4 for a two-level 
atom, with the permittivity and permeability of the plate being again given 
by Eq. (217). It is seen that for a thin plate the barrier is very low. It raises 
with increasing thickness of the plate, reaches a maximal height for some in- 
termediate thickness and then lowers slowly towards the asymptotic half space 
value as the thickness is further increased. 



4.1.3.4 Planar cavity When an atom is placed between two magneto- 
electric plates, then the competing effects of attractive and repulsive interac- 
tion with the two plates can result in the formation of a trapping potential 
[491,497]. Let us model a magneto-electric planar cavity by two identical half 
spaces of permittivity s(u) and permeability fi(u) which are separated by 
a distance d. 29 Substitution of the Fresnel reflection coefficients (B.10) into 

29 Purely electric planar cavities have been modeled with various degrees of detail, 
e.g., by two perfectly conducting plates [227,290,293,298-301,426], two electric half 
spaces [426], or even two electric plates of finite thickness [310,311]. 
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Fig. 5. The potential of a ground-state two-level atom placed between two (a) mag- 
neto-dielectric (all parameters as in Fig. 3), (b) purely dielectric [//(u>) = 1, other 
parameters as in (a)], and (c) purely magnetic [e(uj) = 1, other parameters as in 
(a)] half spaces of separation d = 15c/cjio, Eq. (223), is shown as a function of the 
position of the atom. 

Eq. (203) yields the potential of an atom placed within a cavity bounded by 
the half spaces: 



dq 



q 



e -2bz A _)_ e -2b(d~z A ) 



1 tAjQb-h 



-1 + 2 



q 2 c 2 \ 1 e(ig)6-6i 
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(223) 



As expected, the potential is in general not the sum of the potentials associated 
with the left and right plates separately, as a comparison of Eq. (223) with 
Eq. (213) shows. Clearly, the difference is due to the effect of multiple reflection 
between the two plates which gives rise to the denominators D a , 



1 



e- bd r a+ e- bd 



n=0 



(224) 



recall Eq. (149). 
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The formation of a potential well is illustrated in Fig. 5 where the potentials 
of an atom placed between purely dielectric and purely magnetic plates, are 
also shown. It is seen that the attractive (repulsive) potentials associated with 
each of two purely dielectric (magnetic) plates combine to an infinite potential 
wall (well) at the center of the cavity, while a potential well of finite depth 
can be realized within the cavity in the case of two magneto-dielectric plates 
of sufficiently large permeability. 



4-1-4 Asymptotic power laws 

As we have seen, the dispersive interaction of polarizable/magnetizable objects 
in their ground states can often be given in terms of simple asymptotic power 
laws in the retarded and non-retarded limits. Typical examples are given in 
Tab. 1 where the asymptotic power laws for the dispersion force on an atom 
interacting with another atom [Eqs. (196), (197), (201) and (202)], a small 
sphere [488], thin ring [490], a thin plate [Eqs. (220)-(222)] and a semi-infinite 
half space [Eqs. (214)-(216)], and for the force per unit area between two half 
spaces [372,485], are shown. 

Comparing the dispersion forces between objects of different shapes and sizes, 
it is seen that the signs are always the same, while the leading inverse powers 
are the same or changed by some global power when moving from one row of 
the table to another. This can be understood by assuming that the leading 
inverse power is determined by the contribution to the force which results 
from the two-atom interaction [row (a)] by pairwise summation. Obviously, 
integration of two-atom forces over the (finite) volumes of a small sphere (b) 
or a thin ring (c) does not change the respective power law, while integration 
over an infinite volume lowers the leading inverse power according to the 
number of infinite dimensions. So, the leading inverse powers are lowered by 
two and three for the interaction of an atom with a thin plate of infinite lateral 
extension (d) and a half space (e), respectively. The power laws for the force 
between two half spaces (f) can then be obtained from the atom-half-space 
force (e) by integrating over the three infinite dimensions where integration 
over z lowers the leading inverse powers by one while the trivial integrations 
over x and y yield an infinite force, i.e., a finite force per unit area. It follows 
from the table that many-atom interactions do not change the leading power 
laws resulting from the summation of pairwise interactions, but only modify 
the proportionality factors. 

All dispersion forces in Tab. 1 are seen to be attractive for two polarizable 
objects and repulsive for a polarizable object interacting with a magnetizable 
one. It can further be noted that in the retarded limit the forces decrease more 
rapidly with increasing distance than might be expected from considering only 
the non-retarded limit. Finally, the table shows that the retarded dispersion 
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Table 1 

Asymptotic power laws for the forces between (a) two atoms, (b) an atom and a 
small sphere, (c) an atom and a thin ring, (d) an atom an a thin plate, (e) an atom 
and a half space and (f) for the force per unit area between two half spaces. In the 
table heading, p stands for a polarizable object and m for a magnetizable one. The 
signs + and — denote repulsive and attractive forces, respectively. 



forces between polarizable/polarizable and polarizable/magnetizable objects 
follow the same power laws, while in the non-retarded limit the forces between 
polarizable and magnetizable objects are weaker than those between polariz- 
able objects by two powers in the object separation. This can be understood 
by regarding the forces as being due to the electromagnetic field created by 
the first object interacting with the second. While the electric and magnetic 
far fields created by an oscillating electric dipole display the same distance de- 
pendence, the electric near field (which can interact with a second polarizable 
object) is stronger than the magnetic near field (which interacts with a second 
magnetizable object) by one power in the object separation (giving rise to a 
difference of two powers in second-order perturbation theory). 
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4-2 Excited atoms 



In a first attempt, dispersion forces on atoms in excited energy eigenstates 
can also be derived from perturbative energy shifts. A straightforward gener- 
alization of the calculation outlined in Sec. 4.1.1 to an atom prepared in an 
arbitrary energy eigenstate \m) yields the potential [465,491,492] 



where 



and 



U m (T A ) = U%(r A ) + U* m (r A ) (225) 
£CM = ^ / o °°da 2 Tr[a m (iO-G«(rA,^,iO] (226) 
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[0(z), unit step function] are the off-resonant and resonant contributions to 
the potential, respectively, and 



<y. m (uo) = lim — 
v 1 ^oh^ 



u km -u-ie u km + lo + ie 



(228) 



is the atomic polarizability tensor. Equation (225) [together with Eqs. (226) 
and (227)] obviously reduces to the ground-state result (186) in the special 
case m — 0. Note that the resonant contribution vanishes in the ground state; 
it is only present for an excited atom that can undergo real transitions. Equa- 
tions (225)-(227) which can also be obtained by means of linear-response 
theory [412,413], have been used to calculate the potential of an excited atom 
near a perfectly conducting [412,413], dielectric [413], birefringent dielectric 
half space [428] and an electric cylinder [425]. 

The above mentioned approaches are time-independent and essentially per- 
turbative and inspection of Eqs. (225)-(228) reveals that the application of 
(static) perturbative methods to excited atoms is problematic in several re- 
spects. Firstly, the potential is determined by quantities that are attributed to 
the unperturbed atomic transitions which do not take into account the effect 
of line broadening, whereas in practice finite line widths are observed which 
are known to strongly affect resonant transitions. Secondly, the potential and 
hence also the force remains constant in time; this is not very realistic for 
excited atoms which undergo spontaneous decay with the allowed (dipole-) 
transitions being the same as those entering the potential. And thirdly, per- 
turbation theory does not apply to the case of strong atom-field coupling. 
These problems can be overcome by a dynamical approach to the calculation 
of forces acting on excited atoms. 
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4-2.1 Dynamical approach 



Instead of deriving the dispersion force from an energy shift by some means or 
other, we return to the origin of the force by starting from the Lorentz force 
acting on an atom and calculating its expectation value for a given initial 
state. In particular, when the electromagnetic field is initially in its ground 
state, then this expression yields the sought dispersion force which is genuinely 
time-dependent for atoms initially prepared in excited states. 

Summing the physical momenta m a r a of the particles constituting the atom 
[as given by Eq. (114)], one obtains for the atom as a whole 

mA = PA + /d 3 rP at (r)xB(r). (229) 

Hence, the center-of-mass motion is governed by the Newton equation 

m A r A = F L (230) 
where, according to Eq. (229), the Lorentz force is given by 
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with P a t(r) from Eq. (79). The first term in Eq. (231) can be further evaluated 
by recalling Eq. (107) and using the commutation relations (76). By making 
use of the identity V4p at (r) = — VP at (r) [recall Eq. (79)], one can show that 



1 

2To 



d 3 rp2 t (r),p A 



1 

2To 



d 3 rVp2 t (r) = 



(232) 



and by recalling Eq. (106) and using the definitions (79) and (80), one derives 



E 



a 2m a 



Pc 



d 3 rH a (r)xB(r) \ ,p A 



= V, 



d 3 r 



M at (r) + P at (r)xr A 



B 



(r)}. (233) 



Equations (232) and (233) then imply that Eq. (231) can be written as 



+ 



d_ 
df 



d 3 rP at (r)-E(r) 



d 3 rP at (r)xB(r). 



d 3 r 



M at (r) + P at (r)xr A 



B(r) 



(234) 



It should be mentioned that by using Eqs. (81)-(83) together with the Maxwell 
equations (95) and (97), this equation can be given in the equivalent form 



F L = J d 3 r[p at (r)E(r) +j at (r)xB(r) 



(235) 
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which corresponds to the Eq. (120) used in Sec. 3 as a starting point for 
calculating dispersion forces on bodies. 30 

From Eq. (234), the Lorentz force in long- wavelength approximation can be 
obtained by performing a leading-order expansion in the relative particle co- 
ordinates f a , resulting in 



d-E(r A ) + m-B(r A ) + dxt A -B(r A ] 



+ 



dt 



dxB(r A ) 



(236) 



[recall Eqs. (87) and (88)] where 
d 



dt 



dxB(f/ 



h 



H, d x B(r A )J = d x B(r A ) + d x B(r) 
(h-V A )B(r A ) + B(r A )(y A -i A ) 



r=r A 



(237) 



Discarding all terms proportional to (which are of the order v/c and thus 
negligible for nonrelativistic center-of-mass motion) , as well as the contribution 
from the magnetic interactions, Eq. (236) reduces to 



F L = V 



whose expectation value 



d-E(r) 



d 

+ d^ 



dxB(r) 



(238) 



r=r A 



F = (F L ) (239) 

quite generally provides a basis for calculating electromagnetic forces on non- 
magnetic atoms, including dispersion forces. Needless to say that Eq. (238) is 
valid regardless of the state the atom and the body-assisted field are prepared 
in. 

At this point we recall that, according to Eqs. (44) and (47), the electric 
and the magnetic induction fields are expressed in terms of the dynamical 
variables f\(r, u) and fj(r, u). It is not difficult to prove that in electric-dipole 
approximation, f A (r, uj) obeys the Heisenberg equation of motion 



h 



-kuf A (r, u) + - d- G* (r A , r, u) (240) 



[recall the Hamiltonian (109) together with Eqs. (110), (111) and (119)], whose 
formal solution reads 



f A (r, U,t)= fxfree(r, UJ,t) + f As ourcc(r, U, t) 



(241) 



30 Note that the field created by the atom only gives rise to internal forces, so that 
one may equivalently write Eq. (235) with the total fields £ (r) and B(r) [Eq. (93)] 
instead of E(r) and B(r). 
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where 

f Afrec (r, W ,t)=e" i ^^')f A (r,^) (242) 

and 

fAsourcc(r,^,t) = l - rdre- i ^d(r)-G A [r A (r),r,^] (243) 
a Jt 

(to, initial time), respectively, determine the free-field parts E free (r, cu, t) and 
B frec (r, ou, t) and the source-field parts Es urce( r > u, t) and Bs urce( r , u, t) of the 
electric and the induction field in the ui domain. Substitution of Eqs. (241)- 
(243) together with Eqs. (44) and (47) into Eq. (239) together with Eq. (238) 
and use of Eq. (45) leads to the following expression for the mean force 
[465,492]: 

F(t) = Ffree(*) + 

*- source 

(t) (244) 

with 



p frcc(t) = jf <M V(d(t)-E frcc (r,a;,t) 



d 

+ dt 



d(t)xB iree {r,cu,t)) +C.c. (245) 



r=r A (t) 



and 

F S ourcc(t) — F source (t) + F sou ^ cc (t) (246) 

where the components 



-pel 

source 



(t) = ftfW fdre-M*-) 

[ TV JO Jto 

xV(d(t)-ImG[r,r A (r),^d(r)) + C.c. (247) 



and 



dt Jto 

x (d(t) x (V x Im G[r, t a {t),u]) -d(r)) 1 + C.c. (248) 

J r=f A (t) 

are related to the source-field parts of the electric and the induction field, 
respectively. 

While Eq. (244) together with Eqs. (245)-(248) gives the force on a (non- 
magnetic) atom subject to an arbitrary electromagnetic field, the pure disper- 
sion force can be obtained by considering the case where the (body-assisted) 
electromagnetic field is initially prepared in the ground state |{0}) so that 

({0}|---E frcc (r^,t)|{0}) = ({0}|---B frcc (r,^ ! t)|{0}) = (249) 
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[recall Eq. (63)] which implies that Ff ree (t) = 0. Hence, Eq. (244) simply reduces 
to 

F(t) = F sourcc (t) (250) 

in this case. In particular, for chosen atomic position, may be regarded as 
a time-independent c-number parameter [r A (t) i— > r^], so that the expectation 
values to be taken in Eqs. (247) and (248) only refer to the internal state 
of the atom. It should be pointed out that the concept is not restricted to 
the calculation of the mean force but can be extended to higher-order force 
moments (for a discussion of force fluctuations, see also Ref. [294]). 

4-2.2 Weak atom-field coupling 

The remaining task now consists in the determination of the dipole-dipole 
correlation function 

dmnd m 'n' {^A mn (t) A m i n i (t)^ (251) 

m,nm',ri 

in Eqs. (247) and (248) [A mn = \m)(n\, recall Eq. (111)]. To that end, the 
problem of the internal atomic dynamics must be solved. Let first consider 
the case of weak atom-field coupling where the Markov approximation can 
by used to considerably simplify the problem. Under the assumption that the 
relevant atomic transition frequencies are well separated from one another so 
that diagonal and off-diagonal density matrix elements evolve independently, 
application of the quantum-regression theorem (see, e.g., Ref. [450]) yields the 
familiar result 

(A mn (t)A m , n ,(T)) = 5 nm , (A mn/ (T)) e {^A)-\T m (r A ) + T n{rA)] ,2} ( t-r) ^ 

(t>r, n). Here, 

0J mn (r A ) = uj mn + 6u m (r A ) - 8u n (r A ) (253) 

are the atomic transition frequencies including the position-dependent energy- 
level shifts 31 

<^ m (r A ) = £5u4(r A ), (254) 

k 

5u m {r A ) = —V \ dcucj — 255 

nh Jo oo mk (r A ) - uj 



The Lamb shifts observed in free space are thought of as being already included 
in the frequencies to mn . 
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(cf. also Refs. [493-495]) which are due to the interaction of the atom with 
the body-assisted electromagnetic field, and similarly, 

r m (r A ) = £rt(r A ), (256) 
k 

r m( r A) = ^Q[^ m k(r A )}^ m k{rA)d km dmG[r A ,r A ,uj mk (r A )]-d mk (257) 

are the position-dependent level widths. Note that the position-dependent 
energy shifts h8cu m (r A ) as given by Eq. (254) together with Eq. (255) reduce 
to those obtained by leading-order perturbation theory, Eqs. (225)-(227), if 
the frequency shifts in the denominator on the r.h.s. of Eq. (255) are ignored. 

Substituting Eqs. (251) and (252) into Eqs. (246)-(248), one can then show 
that the force on an atom that is initially prepared in an arbitrary state can 
be represented in the form [465,491,492] 

F(i) = 5> nm (f)F ron (r A ) (258) 

m,n 

where the atomic density matrix elements a nm (t) = (A mn (t)) solve the intra- 
atomic master equation together with the respective initial condition, and we 
have 



F„ m (r A ) = F^(r A ) + F^(r A ) + + KSTM (259) 

with the various electric/magnetic, off-resonant/resonant force components 
being given as follows: 

F mn M = — J o d « 

run -iO]-G«(r A ,r,iO}) r=r , (260) 



F^M=|loE©(Wn fc )J^ m *M 
k 

x {Vd mk -G {1) [r,r A ,Q mnk {r A )]-d kn \ + C.c, (261) 



2ir Jo 



d££ 2 Tr< 



j| «mn( r A,-lO 



X 



VxG (1) (r,r A ,iO 



, (262) 



r=r A 



66 



F^f- r (r A ) = HoJ2®( lI} nk)timn(rA)n mnk (r A ) 

k 

(d mk x{VxG {1) {r,r A ,n mnk (r A )]-d kn }) r=r + C.c. (263) 



k 

X 



Here, Vt mnk {v A ) and ai. mn (r A ,u), respectively, are the complex atomic transi- 
tion frequencies and the generalized polarizability tensor: 



tt m nk(r A ) = u nk {r A ) + \[T m {v A ) + r fe (r A )]/2, (264) 

(265) 



a mn (r A , to) = - ^2 
n k 



d mk d kn d kn d mk 



-VL mnk (T A )-U -Kmk M + U 



Equations (258)-(263) show that the force components F mn (r A ) (m^n) asso- 
ciated with (non- vanishing) off-diagonal elements of the atomic density matrix 
contain contributions arising from the interaction of the atom with both the 
electric and the magnetic field, where the magnetic force components display 
a vector structure which is entirely different from that of the electric ones. 
Since under the assumptions made, diagonal and off-diagonal density matrix 
elements are not coupled to each other so that 

(m ^ n), force components associated with off-diagonal density-matrix ele- 
ments can only be observed if the atom is initially prepared in an at least par- 
tially coherent superposition of energy eigenstates. Accordingly, if the atom is 
initially prepared in an incoherent superposition of energy eigenstates, then 
only force components F mm (r A ) which are associated with diagonal density- 
matrix elements and which are electrical by their nature, 

F mm (r A ) = F%Z(r A ) + F^(r A ) = F£ m (r A ) + F^ m (r A ), (267) 

are observed, with the density matrix elements obeying the balance equations 

(r A )a mm (t) + r™(r A )a kk (t). (268) 
k 

In particular, when the level shifts and broadenings are neglected, then the 
force components F^ m (i\4) and F r mm (T A ) as follow from Eqs. (260) and (261) 
obviously reduce to those that are obtained from the perturbative poten- 
tial (225)-(227) by means of Eq. (183). Note that the gradient in Eqs. (260) 
and (261) acts only on the Green tensor and not on the additional position- 
dependent quantities, so that this result cannot be derived from a potential in 
the usual way. Since the force components associated with excited-state den- 
sity matrix elements are transient, they are only observable on time scales of 
the order of magnitude of the respective decay times T^(r A ) which are known 
to sensitively depend on the atomic position [450]. Needless to say that the 
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Fig. 6. The resonant component of the force on a two-level atom in the upper state 
placed in front of a dielectric half space, Eq. (271), is shown as a function of the 
unperturbed transition frequency wio (solid line) (ujp e /uj^ e = 0.75, ^ e /u)i e = 0.01, 
w Tet d oi + ( d oi-e z ) 2 ]/37rfe c 3 = 10~ 7 , z A /\ Te = 0.0075, A Te = 2Trc/u Te ). For com- 
parison, both the perturbative result as obtained from Eq. (227) (dashed lines) and 
the separate effects of level shifting (dotted lines) and level broadening (dash-dotted 
lines), are shown. 

force F(t) that acts on an initially excited atom approaches the ground-state 
force F 00 (r^) after sufficiently long times, lim t _ i . 00 (F(£)) = F o(i"a). 

In order to illustrate the effect of the body-induced level shifting and broaden- 
ing on the force, let us consider a two-level atom which is situated at distance 
z A very close to a dielectric half space. By means of the respective Green ten- 
sor (App. B), it turns out that in the non-retarded limit the shift and width 
of the transition frequency are determined by 



5uj(z a ) = 8u 1 {z a )-8uq(z a ) 



doi + (d i-e 2 ) 2 \e[u)va + 8u(z A )}\ 2 - 1 



and 



Y{z A ) = T 1 (z A ) = 



d 2 

u oi 



32irheoz A \e\wio + 5uj{z a )} + 1| 



l i-e 2 ) 2 Initio + 5u(z A )\ 



8nheoz A \e[uio + 8uj(z a )] + 1| 



(269) 



(270) 



respectively, where the transition-dipole matrix element has been assumed 
to be real and the (small) off-resonant contribution to the frequency shift 
has been omitted. Note that due to the appearance of the frequency shift 
on the r.h.s. of Eq. (269), this equation determines the shift only implicitly. 
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The dominant contribution to the force on the atom in the upper state is the 
resonant one, i.e., Fu(ta) — F t u (ta)- Substituting the half-space Green tensor 
(App. B) into Eqs. (261) and (267), one can show that [F r n (r A ) = F 1 r 1 (^ A )e 2 ] 
[465,491,492] 

3[d 2 1 + (d 01 .e 2 ) 2 ] |^ 110 (^)]| 2 -1 (97U 
Pll[ZA) ~ 32nhe z\ \e[Sl lw (z A )] + 1| 2 1 ' 

where, according to Eq. (264), 

n 110 (z A ) = u>w(z A ) + iT(z A )/2 = uj 10 + Suj(z a ) + iT(z A )/2. (272) 

In particular, for a medium whose permittivity is of Drude-Lorentz type, 
Eq. (217) leads to (7 e ,T<a;Te) 

u 2 

e[n 110 (z A )] = 1 + — ^y-— * e (273) 

^Te - ^loM - 1)P(ZA) + le\U W (z A ) 

showing that the absorption parameter of the half-space medium, 7 e , is re- 
placed with the total absorption parameter, i.e., the sum of 7 e and the spon- 
taneous-decay constant T(z A ) of the atom. Figure 6 displays the resonant 
component of the force on a two-level atom in the upper state placed near a 
(single-resonance) dielectric half space as a function of the unperturbed tran- 
sition frequency u\q. It is seen that in the vicinity of the (surface-plasmon) 



frequency c^s — y^Te + an enhanced force is observed which is attrac- 

tive (repulsive) for red (blue) detuned atomic transition frequencies uiq < ujg 
(cuio > us) — a result already known from perturbation theory [413]. However, 
it is also seen that due to body-induced level shifting and broadening the ab- 
solute value of the force can be noticeably reduced. Interestingly, the positions 
of the extrema of the force remain nearly unchanged, because level shifting 
and broadening give rise to competing effects that almost cancel. 

The calculation of the off-resonant component of the force, Eq. (267) together 
with Eq. (260), leads to [F^r^) = F°{{z A )e z \ 

3[d 2 1 + (d or e 2 ) 2 ] r e(iO - 1 



32n 2 he z\ Jo s e(i£) + 1 

u w (z A ) ^ (^)+e 2 + r 2 (^)/4 
^i 2 o(^) + [e + r(^)/2] 2 ^ 2 (^) + [e-r(^)/2] 2 - 1 ' 

Equation (274) reveals that the off-resonant component of the force is only 
weakly influenced by the level broadening [the leading-order dependence be- 
ing 0(r 2 )] which is in agreement with the physical requirement that the vir- 
tual emission and absorption processes governing the off-resonant component 
should be only weakly affected by decay-induced broadening. Formally, the 
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1.12 1.13 1.14 

LO\q/ LOTe 

Fig. 7. The off-resonant component of the force on a two-level atom in the upper 
state placed in front of a dielectric half space, Eq. (274), is shown as a function of 
the unperturbed transition frequency (solid line), the parameters being the same as 
in Fig. 6. For comparison, the perturbative result is also shown (dashed lines). The 
inset displays the difference between the force with and without consideration of 
level broadening (solid lines). In addition, the same difference is displayed when the 
level shifts are ignored (dashed lines). 



absence of a linear-order term 0(T) is due to the fact that the atomic polariz- 
ability (265) enters the off-resonant force components (260) only in the com- 
bination [a mm (r J 4, i£) + ct mm (rA, — i£)], which could not have been anticipated 
from the perturbative result (226) [where in fact, 0^(1% i£) and a m (i\4, — i£) 
coincide, recall Eq. (228)]. The effects of level shifting and broadening are il- 
lustrated in Fig. 7. Compared to the perturbative result, the frequency shift 
has the effect of raising or lowering the force for uio < us or cui > us, respec- 
tively, whereas the effect of broadening is not visible in the curves. Only by 
plotting the difference between the results with and without broadening, a 
slight reduction of the force becomes visible in the vicinity of c^s where T is 
largest. Since this behavior is generally typical of off-resonant components, the 
perturbative result may be regarded as a good approximation for the forces 
on ground-state atoms where no resonant components are present. 
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4-2.3 Strong atom-field coupling 



Strong atom-field coupling may occur if an initially excited atom interacts 
resonantly with a sharply peaked (quasi-)mode of a body-assisted field, as 
it is observed in cavity-like systems. In this case, the atom-field dynamics 
can no longer be described within the Markov approximation. Typically, a 
single atomic transition is in resonance with such a mode, so the (resonant 
part of the) dispersion force can be studied, to a good approximation, by 
employing the two-level model in rotating- wave approximation with respect to 
the interaction Hamiltonian (119) (see, e.g., Ref. [450]). To be more specific, 
let us consider the interaction of a two-level atom initially prepared in the 
upper state |1) with the body-assisted electromagnetic field in the ground 
state |{0}) and calculate the electric part of the resonant component of the 
force acting on the atom, i.e., in the Schrodinger picture, 

F(t) ~ Mt)|{v[d.E(r)]} r= J^)>, (275) 

with the state vector \ip(t)) [\i/)(t = t )) = |{0})|1)] being represented in the 
form 

+ E m J f da, ^| d ra ■ Gl (r„ r, u ) ■ |l(r„ .1)10). (276) 

It is normalized to unity provided that 

\Mt)\ 2 + du\Mu,t)\ 2 = l (277) 
Jo 

and 

g 2 (v A ,cu) = ^co 2 d 10 -lmG(r A ,r A ,cu)-d 01 . (278) 
Substituting Eq. (276) into the Schrodinger equation 

with H being given according to Eq. (109) together with Eqs. (110), (111) and 
(119), one obtains the following coupled differential equations for ipi(t) and 

i f°° 
M*) = -tEiM*) +i / diug(r A ,uj)Mu,t), (280) 
n Jo 

ip (u, t) = ~ (E + M^ (w, t) + ig(r A , uty^t). (281) 

Equation (281) together with the initial condition ipo(u,t = t ) = can be for- 
mally integrated in a straightforward way. Inserting the result into Eqs. (276) 
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and (280) then yields 
\^(t)) = ^i(*)|{0»|l> + l - V / d 3 r H duo [' dre-^^+^-^^r) 

x dm • G* x (r A , r, uj) ■ \l(r A , u))\0) (282) 

and 

ly 1 POO ft 

Mt) = -i-±Mt)- dwg 2 (r A ,u) dr e-W^'-^r), (283) 
n Jo Jt 

respectively. Combining Eqs. (275) and (282) and making use of the integral 
relation (45), one finds that 

F (i) = ^° f°<W{V [d 10 -Im G (1) (r, r A ,oj)-d 01 ] } 

x / t dr^(0^i(^)e~ w '^ )(t ~ T) +C.c. (284) 



Since so far nothing has been said about the strength of the atom-field cou- 
pling, Eq. (284) gives the electric part of the resonant component of the dis- 
persion force on a two-level atom which is initially prepared in the upper state 
for arbitrary coupling strengths. Let us now approximate that part of the ex- 
citation spectrum of the body-assisted electromagnetic field which may give 
rise to strong atom-field coupling in the resonant transition by a quasi-mode 
(labeled by v) of Lorentzian shape, 

g 2 (r A , uj) = g 2 (r A , uj,) _ ^ + g' 2 (r A , uj) (285) 

[uj uj v ) -+- 7w ^ 

and assume that the effect of the residual part of the field which is described 
by the term g' 2 (r A , uj) is weakly coupled to the atom, so that it can be treated 
in the Markov approximation. From Eq. (283) it then follows that [498] 

fa(t) = e [-^i/h-iSu>[(r A )-T[(r A )/2](t-to)^ (286) 

where 4>\{t) is the solution to the differential equation 

Mt) + {iAcu(r A ) + [ lu - ri(r A )] /2}0x(t) + ^(r^)^) = (287) 
together with the initial conditions 0i(t = to) — 1, <f>i(t — t ) = 0. Here, 



Auj(r A )=uj u — uj' 10 (r A ) and f^r^) = y2Tcy u g 2 (r A ,uJ u ), respectively, are the 
detuning and the vacuum Rabi frequency and 

fo ; M = fol M + V[A^)]4 7 ;/4 (288) 
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and 

r- lM = r lM - ^ [A ^> + 7;/4 , (289) 

respectively, are the shift and width of the upper level associated with the 
residual part of the field where Sui(t a ) and r^r^), are defined according to 
Eqs. (254)-(257) with the shifted transition frequency being given by 



,32 



u>' 10 (t a )=u 10 + 6u' 1 (t a ) (290) 
in place of Eq. (253). Equation (287) can easily be solved to obtain 

Mt) = c + (r A )e n + (rA ^^ + c_(r A )e n -( r ^('-'°) (291) 

where 

= n^-i ( r A) < 292 » 

and 

n±(r A ) = - IjiA^M+^-riM]^} 



Ti^iAcM + ^-riM]^} -fi|(r A ). (293) 

Combination of Eqs. (284), (286) and (291) then leads to the sought force 
[498]: 

F(f) = — r duu 2 s(r A ,u,t-t ){ J v\d 10 -lmG (1 \r,r A ,Lu)-d 01 }) + C.c. 

IT JO I L J J r= r A 

(294) 

with 
s(r A ,v,t) 

e [-r' 1 (r A )+ni(r A )+n + (r A )]t _ e {i[w' 10 (r A )-u;]-ri(r A )/2+fi:;(r A )}t 

= |c+M|2 c-^oM + iriM/2-ifi+M 

e [-r' 1 (r A )+n;(rA)+n_(r A )]t _ e {i[fi/ 10 (p jl )-a;]-r' 1 (p i4 )/2+n* (r A )}t 

+ c;(r A )c_(r A ) c-^oM + inM/2-i^M 

e [-r' 1 (r A )+H* (r A )+Q+(T A )]t _ e {i[u/ 10 (r A )-u;]-r;(r^)/2+fil(r A )}t 

+ c*_(v A )c + (v A ) c-^oM + inM/2-i^M 

e [-r' 1 (r A )+n*_(r A )+n_(r A )]t _ e {i[d,[ (r A )- W ]-r' 1 (r 4 )/2+fil(r A )}t 

+ |c ~ M|2 c-^oM + irUr^/2-i^M • (295) 

Let us first make contact with result obtained in the limit of weak atom-field 
coupling where the first term under the square root in Eq. (293) is much larger 



32 Note that contrary to Eq. (253), the ground-state shift is absent here as a conse- 
quence of the rotating-wave approximation. 
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than the second one. This is the case when for given transition dipole moment, 
the spectrum of the field in the resonance region is sufficiently flat, 

lv > 2fi R (r A ) (296) 

or when the atomic transition is sufficiently far detuned from the field reso- 
nance, 

|Aw(r A )|»2n|(r A )/ 7l/ . (297) 
By means of Taylor expansion it can then be shown that 



Q±(r A ) ~ < 



-iAu(r A ) - [ 7 „ - r[(r A )}/2 : 

i^|(r A ) Au(r A ) tt R (r A ) lu 



I 4 [A W (r,)F + ^/4 8 [Acu(r A )] 2 + 7 2/4' 

(298) 

so c+(r A ) ~ 1, c_(r A ) ~0 and Eq. (284) [together with Eqs. (286) and (291)] 
approximates to [498] 

1/ roo Vdio-Im G^(r, r A ,u)-d 01 

F(t) = e- r ^><«-*»> ^ / dc , \ K = r - + C.c. 

7T Jo w - ^io(r A ) - iri(r A )/2 

~ e -ri(rA)(t-to) Fi ( rA ) ( 299 ) 

with 

F 1 (r,)=/ 1 ofi? (r,){Vd r G (1) [r,r Al fl 10 M]-d 01 } + C.c. (300) 

and 

n 10 (r A ) = Q 10 (r A ) + ir!(r A )/2 (301) 

which corresponds to the term an(t)Ff{ r (r A ) in Eqs. (258) and (259) with 
Fi^r^) being given according to Eq. (261). 

The strong-coupling limit is realized if the spectrum of the field features a 
sufficiently sharp peak and the atomic transition is near resonant with this 
peak, such that 

7„ <2n R (r A ) and \Au(t a )\ < 2ft| (r A )/ 7l/ . (302) 
In this case, the square root in Eq. (293) becomes approximately real, 

Q ± (r A ) ~ -\ {iAu;(r A ) + \\ lv - T' X (? A )\\ T ^(r A ) (303) 

where 

fiM = v^rM + [A^MP - [7, - rUr A )] 2 /4 (304) 
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so that the coefficients c±(r A ) [Eq. (292)] can be given in the form 



C±M = 20M ■ (305) 

Substituting Eqs. (303) and (305) into Eq. (294) [together with Eq. (295)], 
one finds that for real dipole matrix elements F(t) approximates to [498] 

F(t) = 2 e -^+ r 'i(^)K*-*o)/2 sin 2[ fi(r . A) ( t _ tQ y 2 ] 

ni(r A) + [a,np - [ 7 ,-riM]V4 FlM (306) 

where Fi(r^) is given according to Eq. (300) with 

Q > w (T A )=u > w (r A ) + iT' 1 (T A )/2 (307) 

in place of Eq. (301). Note that [y v + T' 1 (r A )]/2 = 7^/2 if | Au(t a )\ < 7„/2 and 
[7, + ri(r A )]/2 = rx(r A )/2 if 7 „/2 < | Acu(r A ) | < 2ft|(r A )/ 7i ,. 

Comparing Eq. (306) with Eq. (299), we see that while the resonant compo- 
nent of the force in the limit of weak atom-field coupling simply exponentially 
decreases as a function of time, Rabi oscillations of the force are typically 
observed in the strong-coupling limit — in agreement with the well-known fea- 
tures of spontaneous emission in the two coupling regimes. As a consequence 
of the appearance of Rabi oscillations, the magnitude of the force changes 
periodically; for appropriate spatial structure of the resonantly interacting 
quasi-mode, the atom may be trapped with the trap being set by the atom 
itself, cf. also Refs. [499,500]. Rabi oscillations do not occur if the system is 
initially prepared in a dressed state; in this case the (exponentially decaying) 
force is simply given by the gradient of the position-dependent part ±hQ(r A )/2 
[recall Eq. (304)] of the respective dressed-state energy [501,502]. 



5 Concluding remarks 

Dispersion forces are a particular signature of the quantum nature of the inter- 
action of matter with the electromagnetic field. As soon as the interacting mat- 
ter consists of a large number of elementary atomic particles, exact microscopic 
calculations become very involved. Therefore most theoretical approaches to 
dispersion forces make use of assumptions typical of macroscopic electrody- 
namics by introducing — sooner oder later — familiar macroscopic concepts such 
as boundary conditions at surfaces of discontinuity and/or constitutive rela- 
tions averaged over a sufficiently large number of the elementary constituents 
of the respective material objects. Macroscopic electrodynamics whose appli- 
cability surprisingly ranges even to nano-structures, has the benefit of being 
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universally valid, because it uses only very general physical properties, without 
the need of involved ab initio calculations. Moreover, all the relevant quan- 
tities used for characterizing the material objects can easily be inferred from 
measurements. This concept does not only apply to classical electrodynamics 
but also to QED. Macroscopic QED has been well elaborated for the case 
of locally responding media described in terms of complex- valued, position- 
and frequency-dependent permittivities and permeabilities. It can be extended 
to arbitrary linear media, including spatially dispersing media, since the de- 
scription of the quantized field in terms of current densities and the Green 
tensor associated with the macroscopic Maxwell equations is independent of 
the particular medium description. When supplemented with standard atom- 
field coupling terms, a powerful tool for studying medium-assisted quantum 
effects in QED is obtained. Clearly, the applicability of the theory is restricted 
to body-body and body-atom separations that are sufficiently large compared 
with the length scale on which the atomistic structure of the bodies begins to 
play a role. 

In particular, the so established macroscopic QED provides a unified approach 
to the various types of dispersion forces — an approach which incorporates the 
benefits of normal-mode and linear-response approaches, while exactly taking 
into account real material properties. In particular, dispersion forces between 
electrically neutral, unpolarized and unmagnetized ground-state bodies sim- 
ply reflect the forces which are due to the action of the fluctuating body- 
assisted electromagnetic vacuum on the fluctuating charge and current den- 
sities of the bodies. Since all the relevant characteristics of the bodies enter 
the so obtained force formulas via the Green tensor of the associated macro- 
scopic Maxwell equations, they are valid for arbitrary (linear) bodies. Both 
the Casimir stress and the Casimir force density can thus be introduced in 
a natural way. Moreover, by appropriate Born-series expansions of the Green 
tensor, relations between dispersion forces on bodies and dispersion forces on 
atoms can be established which clearly demonstrate the common origin of all 
these forces. In particular, the force on an atom in the presence of arbitrary 
bodies as well as the force between two atoms can be obtained as limiting 
cases of the body-body force. 

In this article we have restricted our attention to ground-state bodies, i.e, 
to bodies that at zero temperature interact with the electromagnetic field 
where the effect of dispersion forces is purely quantum by nature. As outlined, 
an extension of the central results to include equilibrium systems at finite 
temperatures can be obtained in a straightforward way by simply replacing 
the vacuum averages by thermal averages. In this context it should be pointed 
out that the central assumption of linear-response theories according to which 
the thermal average of the field fluctuations is related to the imaginary part 
of the field response function, is explicitly fulfilled within the framework of 
macroscopic QED. 
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As we have seen, dispersion forces on ground-state atoms turn out to be lim- 
iting cases of dispersion forces on macroscopic bodies, so the corresponding 
formulas can be obtained without explicitly addressing the underlying atom- 
field interaction. Of course, they can also be derived by explicitly solving the 
quantum-mechanical problem of individual atoms interacting with the electro- 
magnetic field, with the presence of macroscopic bodies being again described 
within the framework of macroscopic QED. For ground-state atoms where 
only virtual transitions occur, this leads to results that agree with the ones 
obtained from the purely macroscopic approach, as expected. In fact, explic- 
itly addressing the atom-field interaction is more flexible, because it can also 
be applied to non-equilibrium systems, such as initially excited atoms where 
also real transitions are involved in the atom-field interaction. In this case, 
a dynamical description is in general preferred to be employed, leading to 
time-dependent expressions for the forces, according to the temporal evolu- 
tion of the atomic quantum state. In particular for weak atom-field coupling, 
the force on an initially excited atom is a sum of components whose tempo- 
ral evolution follows that of the associated atomic density matrix elements 
which is in turn governed by the familiar master equation of an atomic system 
undergoing radiative damping. For strong atom-field coupling, damped Rabi 
oscillations may occur which periodically change the magnitude of the force. 
The dynamical approach could serve as a starting point for studying disper- 
sion forces on bodies that are not in thermal equilibrium, by appropriately 
modeling such bodies as collections of excited atoms [281,282]. 

Including linearly responding bodies in macroscopic QED has the advantage 
that from the very beginning of all calculations the effect of the bodies is taken 
into account in a consistent manner, without the need to specify the properties 
of the bodies at an early stage. In this way very general results of broad appli- 
cability can be obtained. This naturally applies not only to dispersion forces, 
but also to other quantum phenomena of radiation-matter interaction which 
are strongly influenced by the presence of macroscopic bodies — phenomena 
that may be subsumed under the term Casimir effect in the broadest sense 
of the word. Typical examples are the enhancement and inhibition of sponta- 
neous emission, resonant energy transfer between atoms or molecules and the 
wide field of cavity-QED effects. 
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A Overview over scenarios 



In order to provide for an overview over the various theoretical works on 
dispersion forces, references containing work on body-body, atom-body and 
atom-atom forces are given in separate tables. In the tables, the references are 
structured according to the scenarios considered, regardless of the methods 
used to address these scenarios. 



Material — > 
Geometry [ 


Perfect cond. 


Electric 


Magneto- 
electric 


Half space 
+ half space 


[242,314-316] 
[322,385,386] 
[397,398] 

[391] T 
[318-320] L 


[175,204,205,208,239-241,281,282] 
[323,325,230-232,236,321,326,329] 
[362-365,370,371,374,378,379,382] 
[387,388,399,433,436-438], [324] T 
[328,330-332,359-361,366-369] T 
[389] L , [348,390] Tq 


[333,334] 
[336,372] 
[393,485] 
[335] T 
[337] T 
[339] T 


Half space 
+ plate 






[445,446] 
[487,480] 


Half space 
+ sphere 


[377,397,398] 


[231,436], [330] T 




Half space 
+ cylinder 


[343,377] 






Plate 
+ plate 




[327,373,381,396,435,483,484] 
[347] T i 




Plate 

+ sphere 




[345] 




Sphere 
+ sphere 


[393,397] 


[235,237,344] 





Table A.l 

Schedular summary of references associated with theoretical work on body-body 
forces. The superscripts indicate that finite temperature (T), lateral forces (L) 
and/or torques (Tq) are included. 
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Material — ► 


Perfect conductor 


Electric 


Magneto- 


Geometry J, 


electric 




[107,177,189-191] 


[174,175,189,190,194,196-199] 


[416,443] 




[245,192,247-251] 


[201-203,206-209,211-213,216] 


[491,492] 




[290,294,404,405] 


[217,219,220,252,305,306,405] 


[497], [309] M 


Half space 


[407], [287,295] E 

Tonr 4 on a ai olE 

[296,304,412] 
[413] , [269] 
[284,408] T 


[407,414,415,417,419,433,441] 

r A nnl Tine 010 ooi ooo ori/ilE 

[490], [195, 2 18, 281, 282, 304] 
[30 / ,30o, 413,41o, 42o,43yj 
[465,491,492] E , [406] M , [360] T 

[ool,4Uy ,41U,44zJ 




Plate 




[444] T 


[491,497] 


Sphere 


|221,42UJ 


rooi ooo onct a i /i /ioi /ioo /i/inl 
[zzl— zz3,3yo,414,4zl,4zz,44UJ 


\ A Q Q~\ 

|4oo| 


Cylinder 




[224-226,414,423,424,440] 
|4y3-4yoJ, [42 5J , [444J 




Planar 


[00^7 onn oin /lo^l 
[22 / ,2y(J,310,426J 


[310,426], [311] 


[491,49 (\ 


Cavity 


|Zyi,Zy (— olllj 
izyzj , izyoi 






Spher. cav. 




[228,229,427,440], [312] E 




Cyl. cav. 




[424,440] 




Parab. cav. 


[302,303] 







Table A. 2 

Schedular summary of references associated with theoretical work on atom-body 
forces. Unless otherwise stated, nonmagnetic ground-state atoms are considered. 
The superscripts indicate that excited atoms (E), magnetic atoms (M) and/or finite 
temperature (T) are included. 
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Material — > 

(rpnmpi'T'v 


Perfect conductor 


Electric 


Magneto- 
electric 


Free space 


[107,164,168-170,174,175,181-183,243-257,265,404,433,489,490] 
[276-282] E ,[179,184] E - M , [178,180,258,266,268-272,274,489] M 
[273] M ' T , [283-286,360,361,408,447] T 


Bulk medium 




[429] 


[467,488,503] 
[446] M 


Half Space 


[214,287-289], [430] T 


[215,216,267,415,431] 


[467,489] 


Planar cavity 


[411], [430] T 


[267], [432] T 





Table A. 3 

Schedular summary of references associated with theoretical work on atom-atom 
forces, possibly in the presence of bodies. Unless otherwise stated, nonmagnetic 
ground-state atoms are considered. The superscripts indicate that excited atoms 
(E), magnetic atoms (M) and/or finite temperature (T) are included. 



B Green tensors 



The Green tensor in free space is given by [448] 



1 (c 



a(£p/c)l-b(£p/c)e p e p (B.l) 



(p = r - r'; p = \p\ ; e p = p/p) where 

a(x) = 1 + x + x 2 , b(x) = 3 + 3x + x 2 . 



(B.2) 



The scattering Green tensor for the planar magneto-electric structure charac- 
terized by Eqs. (143) and (144) is given by [482,504] 

G (1) (r, r', iO = J d 2 q e 1 ^^ G (1) (q, z, z', i£) for < z, z' < d (B.3) 

(q _L e z ) with 



G^(q,z,z',iO 



8tt 2 6 J^J D, 



~2bd 



cr=s,p i 



e+e+e- 6 ^- 2 ') + e-e-e b(z ~ z,) 



1 



■ (B.4) 
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Here, b and D a are defined by Eqs. (148) and (149), respectively, 

ef = e q xe z , e± = -^-(iqe z ± be q ) (B.5) 
(e g = q/g, g = |q|) with 

k = £^r(iOMiO (B.6) 

are the polarization vectors for s- and p-polarized waves propagating in the 
positive (+) and negative (— ) z-directions, and = r^(£,q) with a — s,p 
describe the reflection of these waves at the right (+) and left (— ) walls, 
respectively. 

In particular, assume that both walls are multi-slab magneto-electrics con- 
sisting of N± homogeneous layers of thicknesses d?± (j = 1, . . . , N±) with d^* 
= oo, permittivity e±(oj) and permeability n±(u). In this case the reflection 
coefficients can be obtained from the recurrence relations [r a ± =r° ± ; d = d°± ; 
e(u)= 4M; A*H= A*±H] 

, _ (^ + v ± - ^ +1 ) + 04 +1 % + iMh g^j^rg fB7 x 

rJ (4 + V ± - Mil + + e ±V± l ) e-^^rg 1 , . 

rp± (4 + V ± + 4&i +1 ) + (4 + V ± - 4&i +1 ) e-^^Vg 1 1 ' j 



(j = 0, . . . , iV± — 1) with = where 



*4 = y| 400/400 + <z 2 - (b.9) 

For single, semi-infinite slabs, Eqs. (B.7) and (B.8) reduce to the Fresnel co- 
efficients 

u\b — ub\ e\b — eb\ 

r -=tCT- r ^ = t^i- (R10) 
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